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ABSTRACT 


Five base-line calibrations of the Motorola Mini-Ranger 
III (MRS III) shortwave positioning system were performed 
across Puget Sound, Washington, in January 1984. Two MRS 
III codes were calibrated over measured distances of 1061.2, 
2417.5, 4083.1, 7016.8, and 9861.1 meters. Reference lengths 
were obtained to submeter accuracy using a Tellurometer, 
Model CA-1000. 

Results of the data analysis suggest that an optimum 
base-line length of about 2,500 meters provides a character- 
istic calibration curve that adequately depicts all equip- 
ment configurations for a typical MRS III. The calibration 
curve was best characterized by high-gain antenna data for 
two Mini-Ranger codes. The high-gain antenna also exhibited 
the least sample deviation of the three antenna types cali- 
brated. Range average was found to be insignificant.  Re- 
gression curves applied to the data showed no change in 
quality of fit above a power of three and variances of 0.5 


meters indicateda fair to good fit of the data. 


TA 


IT: 


TABLETOF CONTENTS 


INDRODUCTION ο ο ο WES... n ME ΝΥ ΠΡ’. 
А. BACKGROUND TINEORMATION 1 2.2. (ELE «ον ο 
1. SHEj Positioning Systems 0 .0 scopo MSS BD 
2. Theory ef Seperation... ae MES dul 
3. Field Uses of SHF Systems ..... νυ eae ct M r2 
Ap “ΜΟΙ OoOons s Be ТЕН! 
В. SYSsRRmMUERPDORS ο S Vu EE E ome e 
Ἱ. SystematiesRange Errors ..... 9 v 94 ь г ΜΜΑ 
2. Signad Sëremetbku asi an Imdieater ..........19 
C. NOS CALZITBRATIONWBRGSGDDURESU. IIR a TA 
1. Ваве-ІБіпе Саймокесіопе .......... те... m 21 
2.  Crdtbudeal Systemoseheoks- ss. Сар 
3. Noncrdstical “Seis tema Песке s - 26 
Ds  POCUS ОББ БББ ЕНК ο s ee КИЕ Ve. MIO 
l.  TheaEuture Or SNP ES. IMEEM όνος ГА 
2. Improving MRS III Calibration Methods ....29 
3. Present Calibration Problems ............ odi 
MINI-RANGER III POSTTIONING SYSEEMSE SS И oe 
A. SYSTEM COMPONE oa νο πο... 33 
l. Receiver-Transmitter Assembly ............ 23 
2. Shore Reference Stations ..... ο ο ο. 35 
3. Antennae man ET КЕГІ” ος 36 
4. Range Comnsole 22222222204 а Ы 41 
B. -MRS III SYSTEM OPTIONS TE O P S P D ΙΡ 43 
l. Range Average 下 к BE 43 


2. Signal SL ЕШШ От ................. 44 


ID E SEISRETMENITSDESPIGNOSANDODATA REDUCTION .............. 47 
J OP p Vivo’ . е, „л a s s e С PEERS 47 
Ελ απο ο. ο... E. NE esa 47 
CA ACOUES LEIPTON AND PROCESSING ο νο Ç s < v e ç e e e ea 48 
Jo U TOT roOnme t eri MeSSUKremeonmtS ,................ 48 
Ze Mint Ranger Data AGquUISTUYOD νεος ο ο ον ο ο ον Б 
3. Mini-Ranger Data Processing ............... p 
J υπ οπα παμε en Me ο «ον νο ο σος 59 
ER 和 RON ND CONCEUSTONS 2..2.................... 80 
"ΕΗ '''. ο δε е ои 25.54. V XM Comememene 58 
Nee ....rs see es es». ος ος аб 
επι капы Е та РТО БИТИ Оаа. лањ 58 

3. Range Error Versus Base-Line 
Di sl JI t So οὐ. «ος «σος «οσον ς 59 
A Ciyvespitting Besults .......... к 60 
Bem CONCLUS TONS ο... . ο. oe a de i r 64 
ГООО ЕАО ИСО ЗОВ араа... 66 
APPENDIX A: MINI-RANGER III SYSTEM SPECIFICATIONS ...... 517 
ΙΙΙ: RANGE ERROR AND STATISTICAL PLOTS .......... Τι 
ΙΗ TFELEUROMETER DATA ........-................. 114 
КИЛЕШЕ COMPUTER PROGRAMS .............„........... T20 
DD ΕΡΕ Ευ ον ος ος ος ον ος ος όν ο ο ον ὁ 147 
有 148 


TI: 


ΤΙ. 


TY: 


MIR 


Vae 


LIST OF TABLES 


ATTENTUATION ESTIMATES FOR A 3-CM WAVE LENGTH ..... 16 
NOS CALIBRATION REJECTION LIMITS ............. την. 7: 
MINI-RANGER COMPONENTS AND SERIAL NUMBERS ........ . 33 
BASE-LINE ENDPOINTS AND RESPECTIVE LENGTHS ........ 50 
REGRESSION CORDETGPENPSEC UE T 2.2222” 60 
COEFFICIENTS FOR CUBIC REGRESSION CURVES .......... 61 
COMPUTED VARIANCES FOR HIGHER DEGREE CURVES ....... 65 


IUSTO SOIN RSS 


ШУО БШ ЫЛОО  Басе-Бл1пе“ Сайт ргакгтоп Сит уе .....---.e.ecv* 18 
pPuacesshapewDegsadation' and" Rdnge Errors ον ος ος ο σσ 20 
"ΡΕ, “snore "Seat rons, and Range CONSOLE 77... 2777 34 
WE ΕΠ CCor ΑΠΕΕΠΠΗΓ.....«.-»ε«.«νν»5οι.:. oy 
INN σα Τη ΟΞ Sector Antenna .SO eee ee ES 
МКО eect Lon т с=п Sector Antennd νο ος ο οσο σσ ο 5 579 
Medium-Gain Antenna Radiation Patterns ............. 40 
"uwHI-Ranger III Range Console Assembly ............. 42 
ο. ScMEBENNEEHONeHEOF .....e«.. ert etm 45 
"ΕΤΕ" and Statrom Босастопв  ¿ Ç s s < ο ο οσο ος ος ο 49 
Range Error versus Base-Line Distance, Code O ...... 54 
Range Error versus Base-Line Distance, Code l ...... 29 
Cc c SONST arves LOr Codon O m ue. s.s s a s a sa a e e n 62 
οι. ο πο ΙΕ s f Oro Е 007 a Ки... 63 


ACKNOWLEDGEMENTS 


The author thanks Ms. Pat Nyland, Mr. Denny Sutton, Mr. 
Rick Hastings, and the command and crew of NOAA Ship RAINIER. 
In addition, Mr. Russ Eastman of the Electronics Branch, 
Pacific Marine Center, provided valuable assistance during 
data acquisition. Capt. Glen Schaefer, thesis advisor, was 
a source of encouragement and assisted in editing. Lt. Cdr. 
Alan Anderson reviewed the thesis and offered many helpful 
suggestions. Photographs of the MRS III were provided by 
the Motorola Corporation. Mr. Brian P. Fields assisted by 


drafting figures. 


L ONBBEODUCTION 


A. BACKGROUND INFORMATION 

Electronic positioning of hydrographic survey vessels is 
a common practice in most automated UEM chic agencies. 
The National Ocean Service (NOS) is such an organization and 
is responsible for charting the United States territorial 
waters. The liability that accompanies this responsibility 
has prompted NOS to analyze the quality of its positioning 
ου pecifically NOS weuld like. to apply functionally- 
derived range error correctors to range rates obtained from 
superhigh frequency (SHF) positioning systems in real time. 
SHF positioning systems, such as Motorola's Mini-Ranger III 
(MRS III), are used for inshore and harbor surveys of scales 
cC ος than  Portability, accuracy, and automation 
are some of the reasons hydrographers prefer them. 

The systems generally read out directly, in meters, the 
distance or range of the vessel. However, these range rates 
are not error free due to weather changes, attenuation dis- 
tances, changing components in a given configuration, and 
other variables. The errors are determined through base- 
line calibrations where observed ranges are compared against 
a known distance and error curves are generated that display 
range error as a function of signal strength. Range correc- 
tors are then extracted from these curves and applied to the 
received range rates in near-real time. 
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In January 1984, calibrations were performed over five 
measured base lines that extended across Puget Sound, 
Washington. Three different antenna-types and three differ- 
ent Range Average settings were sampled for each base line. 
The results were used to explain the observed relationship 
between range error and signal sere over a Variety of 
conditions and system configurations. 

ITCOUSHPSESSItI)MEHUNSYSLUONS 

SHF positioning systems have played a major role in 
the development of automated hydrography in NOS. In general, 
electronic methods have replaced antiquated manual methods 
for positioning on scales as large as 1:10,000. These micro- 
wave "black boxes" were for years rejected by stubborn "moss- 
back" hydrographers who refused to accept the accuracy 
claims. 

Today's electronic systems have freed the hydro- 
grapher of the tedium of manually acquiring and processing 
data so the hydrographer may concentrate on planning, error 
analysis, data interpretation, and other more important 
functions. MRS III has been a mainstay in west coast sur- 
vey operations in NOS since the early 1970's. 

Though most of the problems encountered with MRS III 
have been solved in the course of the agency's experience 
with the system, the problem of systematic range errors and 
their necessary correction has not been fully resolved. 


Additionally, a new automated hydrographic data acquisition 
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system called SDS III is slated to replace the present HYDRO- 
PLOT/HYDROLOG system in 1987. The SDS III development group 
is concerned about present SHF calibration methods and would 
like to devise a procedure for accurately determining range 
correctors over a variety of conditions and applying them in 
near-real time. 

2. Theory of Operation 

The theory of SHF positioning systems is an exchange 
of pulsed signals, where radar-like signals are transmitted 
from the vessel to an active shore station. Upon interroga- 
tion, the shore station returns a signal to the vessel. 
Distance is directly related to the round-trip travel time 
of the signal which is based on a fixed speed of propagation 
of the electromagnetic waves. Atmospheric effects will vary 
this speed of propagation but corrections can be applied 
given the proper meteorological data. [Ref. l:pp. 45-52] 

In the specialized case of the MRS III, the system 
uses one or more reference stations placed over points of 
known position on shore. Each reference station can trans- 
mit different pulse forms (called codes). These codes can 
be set in the reference stations and then selected for in- 
terrogation by setting a dial on the range console. The 
reference stations are actively interrogated by a receiver- 
transmitter assembly (RT unit) located aboard the vessel. 
The range console (also aboard the vessel) contains a stable 


counter clock that measures the travel time and converts it 


ПЕП 


to a distance. The range console also performs some internal 
smoothing and screening functions on the received signals. 
The resultant range rates are displayed at the front of the 
console. The system comes with three different types of an- 
E cete whose purposes are discussed in greater detail in 
Chapter II. The MRS III operates on a frequency of 5.4 to 
5.6 GHz (versus 9 GHz for most SHF systems). This lower fre- 
quency helps reduce the possibility of interference from 
other radars in the working area. The system has been used 
by NOS for hydrographic surveys ranging from scales of 
1:20,000 to as large as 1:5,000, when used in a hybrid mode 
such as range-azimuth. 

3. Field Uses of SHE Systems 

NOS has been using MRS III's for over a decade in a 
variety of survey applications. Most positioning schemes 
use the MRS III in either a range-range or range-azimuth 
geometry. Field edit positions, velocity casts, and bottom 
samples may also be positioned this way. 

SHF systems are very portable and are well-suited 
for inshore hydrography and shoreline development requiring 
a shallow-draft vessel. Also, the positions of control 
stations in a project area can be verified by performing a 
system check calibration of the MRS III in the field. Sys- 
tem checks are discussed in greater detail later in this 


chapter. 
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The systems are also used in NOS to check lane counts 
of medium-frequency positioning systems such as RAYDIST. In 
NOS, MRS III range rates are recorded at periodic intervals 
to explain or pinpoint where RAYDIST lane jumps (lane losses 
or gains) have occurred. Once lane jumps have been deter- 
mined, whole lane correctors can be added to existing partial 
correctors of the RAYDIST to allow the hydrographer to con- 
tinue survey operations. This procedure has saved consider- 
able amounts of hydrographic data by relieving the 
hydrographer of the need to return to a calibration or ini- 
tialization area. 

SHF systems are practical for positioning a vessel 
using range-azimuth geometry. This geometry uses a shore- 
based observer with a theodolite positioned next to the SHF 
shore station. The observer radios observed azimuths to the 
survey launch. This setup is less probable with a medium- 
frequency system like RAYDIST because the observer's pres- 
ence in the ground plane of the antenna affects the outgoing 
signal significantly. The popularity of range-azimuth geom- 
etry has prompted the recent development of automatic and 
semi-automatic total stations that consist of an SHF unit 
and electronic theodolite. Accuracy is increased by using 
an azimuth vice a second range rate in the position 


computation. 
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4- MRS TII Options 


There are several options that can greatly enhance 
the operation of the MRS III and help reduce the effect of 
systematic range errors: 


a. А signal strength indicator displays the strength of 
each incoming signal. Signal strengths can be recorded 
in real time and used to identify position blunders or 
busts. They can also be used to determine range error 
correctors in near real time. 


b. A range averaging option enables the user to select 
the number of range rates to be averaged before dis- 
playing. This option was originally installed to 
provide more accuracy by increased averaging of the 
data. 


c. A multi-user capability offers 16 different codes. 
The extra codes give the hydrographer flexibility to 
work in areas where other Mini-Rangers are operating. 
Though not directly related to range errors, anomalous 
range rates can sometimes be attributed to other MRS 
III systems located in the survey. 


d. A voltage cut-off switch powers the unit down when a 
minimum of 22 volts is not exceeded. It was found 
that voltages below 22 volts caused significant range 
errors. 

e. The Space Diversity option solves the major problem 
of range holes or null zones by using two RT units. 
The nulling phenomenon is prevented by mounting two 
RT unit antennae at different heights on the survey 


vessel with a switch to select between the antennae, 
until an acceptable signal is received. 


BS SYSTEM ERRORS 
1. Systematic Range Errors 
MRS III, though an accurate system (Appendix A) to 
use, still contains "some idiosyncracies that have remained 
largely unknown and which can lead to serious range error 


problems." [Ее?. 2:р. 5] This statement was made іп lovee 
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Today many of the error sources have been resolved but others 
are as mysterious as ever. Operation of the system can be 
affected by range, meteorological conditions, range holes or 
null zones, antenna types, range averaging, power supply 
voltage, warm-up time, and even the RT unit cable length. 
Only errors due to range, range averaging, and antenna types 
were considered in evaluating the January 1984 calibration 
data. 

Changes in weather within a survey project area can 
cause drift in the MRS III. Temperature changes directly 
affect a unit's circuitry--specifically its magnetron. 

"Very few weather phenomena directly or significantly affect 
‘the range measurement of the Mini-Ranger system because of 
the high frequency of operation--5 GHz. However, meteorolo- 
gical conditions do contribute to variations of the effec- 
tive path length (apparent range) and, more severely, 
atmospheric absorption (attenuation)." [Ref. 3:Section 3.1.2] 

The actual propagation speed of the radio frequency 
(RF) waves is a function of wavelength, air temperature, 
barometric pressure, and vapor pressure. The counter clock 
inside the range console uses a fixed, average value for the 
propagation speed. Deviations from the fixed value are 
determined during a base-line calibration performed in the 
field by the user. 

Phenomena of "high particle density such as rain, 


mist, fog, or photochemical smog" are also responsible for 


ls 


attenuation. A formula (Eqn. 1.1) exists for calculating өте 
attenuation due to rain [Ref. 4:p. ll] where A is attenuation 
in db, r is rainfall rate in mm/hour, and a & b are tabulated 


frequency functions. 
A = ar ΓΤ; 


Signal attenuation reduces the range of operation and 
lowers signal strengths, so calibrations should be performed 
as near to the working area as possible to minimize errors 
due to weather. Thus, any range errors found during the 
base-line calibration will be a function of weather condi- 
tions in the project area. Table I shows attenuation esti- 


mates derived by using Equation l.l. 


TABLE I 


ATTENUATION ESTIMATES FOR A 3-CM WAVE LENGTH 


Precipitation Attenuation 

Rate (mm/hr) ( db/km) 
6.0 0.198 
15-2 07613 
πο. 7 0.726 
226 0.801 
34.3 1: 250 
43.1 1.640 


Systematic range errors are the focus for this thesis. 
Their presence is attributed to the fundamental radar equa- 
tion (Eqn. 1.2), where P' is the intensity of radiated power, 
P is the peak power, G is the antenna gain, and R is the 


distance from the source. 


l6 


P' = PG/R (1.2) 


The effects of the radar equation are shown by plotting range 
error versus distance in kilometers. The same curve form is 
also found in plots of range error versus signal strength. 
In the latter case, the errors result from artificially- 
attenuated signals that simulate greater distances. Thus, 
NOS base-line calibrations of the MRS III are related to the 
radar equation. In a typical NOS calibration curve (Fig. 
1.1), the range error curve becomes exponential below a cer- 
tain critical signal strength threshold. At high signal 
strengths the curve exhibits a near linear form. Casey re- 
fers to the threshold value separating the two curve types 
as the Critical Strength Threshold (CST). NOS calls the CST 
the Minimum Acceptable Signal Strength (MASS) and NOS re- 
jects data below this value. Casey has found that above the 
CST "the errors are distributed with a Gaussian or normal 
Ge On; below the CST the distribution is log normal." 
eer. 2p. 5 

The important point is if the curve were considered 
"well-defined" below the MASS, data acquired with signal 
strengths below this threshold might not be rejected. This 
amounts to a large savings in time and money when viewed 
from a (NOS) fleet aspect. А well-defined curve is one 


whose data exhibit low standard deviations. 
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Range errors are due principally to pulse shape deg- 
radation of the received signals. То determine the amount 
of degradation and the effect on the ‚сна range rate, a 
signal strength integer value is used that defines, through 
calibration, what range corrector (the algebraic inverse of 
range error) should be applied to the range rate. 

2: Signal Strength as an Indicator 

NOS MRS III's are equipped with the optional signal 
strength indicator that displays a "relative number output 
that is directly proportional to the amplitude of the re- 
ceived signal from the reference station."  [Ref. 3:Section 
5.2.3] Signal strengths are measured because tests have 
shown that range error is dependent on signal strength in a 
repeatable manner. Signal strengths are significant.because 
the counter clock, which measures the elapsed time, must 
recognize the leading edge of the received pulse because 
". . . the leading edge shape of the incoming pulse varies 
with the amplitude so that a low amplitude (read signal 
strength) leads to a systematic timing (read range) error." 
(Ref. 4:p. 9] This relationship between elapsed counter time 
and amplitude is best shown by diagram (Fig. 1.2). 

The factors affecting MRS III signal strength that 
are discussed in this thesis are range averaging, base-line 


lengths, and three types of reference station antennae. 
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C. NOS CALIBRATION PROCEDURES 
1. Base-Line Calibrations 
Base-line calibrations are required by NOS (1) at 
the beginning and end of a survey project, (2) at intervals 
not exceeding 2 months during long projects, (3) whenever 
an RF component is changed in either an RT unit or tran- 
sponder, and (4) whenever the difference between prior base- 
line correctors and daily system check correctors exceeds 
specified rejection limits. [Ref. 5:2р. М-2] 
A description of the base-line calibration procedures 
is taken from the PMC OPORDER, Appendix M: 
The base-line calibration is performed by introducing 
signal attenuation with a variable attenuator over a 
known distance. The observed range readings are compared 
to the known distance at various signal strengths to 
determine range errors. A graph is constructed 
which typically shows a nearly-linear portion followed by 
a sharp rise at low signal strength values. A mean error 
is determined for signal strength values corresponding 
to the nearly-linear portion of the graph and its al- 
gebraic inverse is the corrector to be applied to the 
observed range readings. The signal strength value 
corresponding to the sharp upswing portion of the curve 
is the minimum value allowed for hydrographic survey 
Operations. [Ref. 5:pp. M-2 & M-3] 
For a base-line calibration, ten range rates are sampled for 
each signal strength to obtain the mean range. Base-line 
lengths are measured prior to calibration to submeter accu- 
racy using an NOS-approved laser or microwave electronic 
distance measuring instrument. 
NOS presently requires the length of the base line 
to be between 700 and 3,000 meters. The optimum length has 


yet to be determined and is discussed in the thesis. 
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Though the correctors determined from a base-line 
calibration are most often those applied to the final posi- 
tion data, NOS does — that a closing or bracketing 
calibration be performed at the a оа of a survey. 
Additionally, daily critical or noncritical system checks 
are made to monitor any drift in the Mini-Rangers. These 
system checks test for instrument drift and verify or dis- 
pute the most recent base-line correctors. 

2. Critical SS cms nechs 

A critical system check is defined as one in which 
an independent method of position determination is employed 
to check range errors. Noncritical system checks are much 
less time consuming and are performed in conjunction with a 
recent critical system check. This approach was devised by 
NOS to minimize the calibration time and maximize the time 
spent acquiring hydrographic data. Though noncritical sys- 
tem checks require less time, they do not provide an inde- 
pendent check fix to verify the first fix. Whenever the 
difference between a daily system check and base-line cali- 
bration corrector exceeds the rejection limits (Table II) 
and cannot be resolved after a careful check is made, a new 
base-line calibration is required to continue hydrographic 
survey operations. [Ref. 5:pp. M-5 & M-6] 

The following are methods currently approved for 


Critical system checks. 
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NOS CALITBRKRATIONTREJECPION LIMITS 


Survey Scale Rejection Limit (m) 
оО 3 
I: 1079090 5 
1: 20,000 and smaller 10 


a.  Fixed-Point Method 
The fixed-point method is the critical system 
check method preferred by NOS. It involves bringing the 
launch alongside a geodetically-positioned pile, wharf, rock, 
or other point feature and comparing the observed range rates 
with precomputed range rates. It is operationally and com- 
putationally simplest of all the methods and is not restrict- 
ed by visibility in the work area. However, the method is 
not always Өзенде as frequently a point feature is not 
available in the work area. Additionally, this method is 
restricted to launches since deep-drafted survey ships can- 
not afford the grounding risk of the shoal waters surround- 
ing such features. 
p. Navigational Range and Cut-Off Angle Method 
The navigational range and a cut-off angle 
method can be employed when as few as two geodetic control 
stations are in the area. To construct a navigational 
range, one geodetic control station is used as a front (or 
rear) range marker. A second station, or erected signal, 


is used as the other range marker. 
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While the survey vessel is moving along the es- 
tablished navigational range, the angle measured by sextant 
between the front (or rear) range marker and the second geo- 
detic control station is the cut-off angle that marks the 
vessel's position when the cut-off angle agrees with any of 
a number of preselected angles. The preselected angles and 
known navigational range azimuth facilitate the precomputa- 
tion of vessel positions. Knowing the vessel position and 
the geodetic position of the MRS III shore stations, pre- 
dicted range rates are computed for each cut-off angle for 
comparison with the observed range rates. 

c.  Range-Azimuth Method 

The range-azimuth method can produce rapid re- 
sults when a total station (a theodolite and electronic dis- 
tance measuring instrument) is used to acquire the azimuth 
and range to the launch from a geodetic control station. 
Given the necessary onboard computer program, a single shore 
observer radios azimuth and range data to the launch. The 
computed position is converted to a set of predetermined 
MRS III range rates that are compared with the observed 
rates. This method is the most accurate available and 
quickest after the preferred fixed-point method. 

d. Intersection of Two Navigational Ranges Method 

The intersection of two navigational ranges 
method can be accomplished with a minimum of four control 


points. A navigation range is constructed as described in 
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the above method (b). Though this may require considerable 
set-up time, the method is quick and convenient once the 
ranges are in place. To calibrate, the launch steers one 
range and crosses the other. The observed range rates are 
compared with precomputed rates by subtraction. 
e. Three-Point Sextant Resection Method 

The three-point sextant resection method has 
been used extensively in NOS pending the development of more 
efficient methods. It requires a minimum of four visual 
targets and three sextant observers. Simultaneous to the 
reading of the MRS III range rates, horizontal sextant an- 
gles are observed to the four signals. The observed sextant 
angles are input into a calibration program. Though compu- 
tations are performed by an onboard computer, they are still 
tedious and more involved than the simpler subtractions used 
in some other methods. 

f.  Theodolite Intersection 

The theodolite intersection method requires 
three observers, each with a theodolite to turn angles from 
an initial station to the survey launch. The recorded azi- 
muths define a fix and check fix that are input into a cali- 
bration program. Although it is one of the more accurate 
methods, it is rarely used by NOS personnel due to the 
equipment and labor demands. The ensuring computations are 


only slightly faster than for the resection method. 
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3. NoncriticaliSvstoR Checks 


Noncritical system checks are generally quicker than 
critical system checks but do not include an independent 
position determination. Additionally, they cannot check geo- 
detic control in the working area. NOS has approved their 
use because the MRS III has been used in west coast survey 
operations for over a decade and has developed into a relia- 
ple system that usually does not develop flaws which cannot 
pe detected using noncritical system checks in conjunction 
with periodic critical system checks. This fact, coupled 
with the policy of using base-line correctors in lieu of 
separate correctors determined daily, has been the driving 
force behind reducing the requirements for twice-daily criti- 
cal system checks. (Кеке Брукс? 

a. Three-Range Method 

The three-range method checks the operation of 
three separate shore stations by comparing the three possible 
fixes generated by the three range rates. A major constraint 
is that the vessel be positioned to obtain a near 60-degree 
angle of intersection for each of the three lines of position 
(LOPs) to obtain consistently strong fixes for each pair. 
Unlike the other two noncritical system check methods that 
utilize ready-made calibration forms, the three-range method 
requires separate fix computations and an onboard computer. 
The key to using this method 15 realizing het it 15 sensi 


tive to the intersection angles of the LOPs. 
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b. Base-Line Crossing Method 

The base-line crossing method compares a sum of 
two observed range rates to the known base-line length con- 
necting two transponder stations. “Thus, the launch crosses 
between two shore stations and a mark is called when a shore 
observer radios that the crossing launch is on range with the 
two shore stations. To stabilize the update rate of the MRS 
III, the launch steers a range arc of one of the stations. 


The ensuing computations are subtractive and simple to 


perform. 
c. Launch to Launch Method 
The launch to launch method is very quick because 
many launches can be checked at once. Using a critically- 


checked launch, designated the guide launch, other launches 
come alongside with their bows pointed at о ло ида: Lo 
be checked. They do not have to be tied but they should be 
aligned along their beams.  MRS III rates are simultaneously 
recorded and the guide launch correctors applied.  Non-guide 
launches then apply their respective base-line correctors 
and compare to the rejection limits table (Table II). А11 
system check methods, both critical and noncritical, require 
the hydrographer to first apply the most recent base-line 


correctors to the observed rates before comparing them 


against their known values. 
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D. ЕОСОЅ ОЕ ЕБЕЕ Е ЕП 
l1. The Future or μμ μυς 

SHF positioning systems have been available for over 
a decade. They may eventually be replaced by futuristic 
satellite systems such as the NAVSTAR Global Positioning 
System (GPS). The need to locate transponders on shore for 
launch operations could be eliminated along with the some- 
times complex geodetic networks that support a survey. With 
only 6 of the scheduled 18 satellites in orbit, testing 
has already begun to determine GPS accuracies. The entire 
System is expected to be operational in 1987. Other advan- 
tages to this powerful system are its continuous passive 
output, multiple user capabilities, and world-wide coverage. 
It will also free the hydrographer from any constraints due 
to fix geometry or line-of-sight. scum have been con- 
ducted by the Naval Oceanographic Office (МАУОСЕАМО) іп Вау 
St. Louis, Mississippi, to determine GPS's suitability for 
dynamic positioning of a survey launch. Thus far, GPS accu- 
racies of 10 meters have been observed in a dynamic mode. 
[Ref. 6*p.944 

A cost study recently done by NAVOCEANO [Ref. 7:p. 
84] states a savings of nearly four million dollars a year 
for two Navy hydrographic survey ships would be a result of 
using GPS. It appears likely that the GPS system will make 
dinosaurs of not only SHF positioning systems, but medium- 


frequency systems as well. If GPS receivers were made 
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affordable even to fishermen, it would replace LORAN-C and 
Omega. However, there are several reasons why the present 
systems will still be used: 


GPS accuracies for hydrography (a dynamic mode) are not 
yet verified. This fact alone suggests that SHF systems 
will continue to be used in NOS for large-scale 
hydrography. 


GPS may not be able to keep pace with the short fix in- 
terval used to obtain hydrographic positions on large- 
scale surveys. Also, additional computer storage may 
be necessary to accomodate the large amounts of GPS 
data. 


Access to the precise (P) code may be restricted to 
military operations. Since NOS does not directly sup- 
port such operations, it may not be granted access to 


this code. Access to the P code also generates con- 
siderable security problems regarding its use by NOS 
personnel. If access to the P code is limited there 


are other techniques, such as interferometry and 
differential GPS, that will improve the accuracy of 
data obtained with the coarse acquisition (CA) code. 
The initial cost for a GPS receiver is prohibitive to 
the average mariner and may also be too expensive for 
NOS. "Initial estimates of costs are $50,000 for 
each system (excluding the cost of the receivers) 

and $15,000 to $30,000 one-time cost for tailored 
software." [Ref. 7:p. 83] Using this estimated, one 
can imagine the initial cost to NOS for nearly 50 
operating systems. 


СЕКЕ ОУ ЛП ος πι Calibration Methods 
The above statements suggest that MRS III and similar 
systems will still be used during the evolution and testing 
of SDS III. Base-line calibrations, to determine systematic 
range errors, will still be necessary. The present base-line 
calibration methods can be automated and enhanced to acquire 
more data in less time without manual manipulation by the ob- 


server. One such practice is to record calibration data 
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directly into a portable computer. The samegiprojgr man c ш 
acquire the data, for example LOTUS 1-2-3, can be used to 
automatically average the correctors and graph them.  NOAA 
Ship DAVIDSON is presently using this procedure for base- 
line calibrations of its Mini-Ranger Falcon 484, the MRS III 
successor. Possible refinements to the method include a 
mathematical, nonsubjective, determination of the MASS and 
applying curve-fitting routines to obtainya functi Бгн ог 
of the calibration curve. 

A further enhancement would be to tabulate the cor- 
rectors into a corrector file such as the Sound Velocity, 
Tide, Zone, Parameter, and Station Tables that exist pres- 
ently in the Bathymetric Swath Survey System (BSSS). 
Another possibility is to incorporate the functional form 
(equation) of the calibration curve directly into the 
start-up routine of the SURVEY or HYDROPLOT acquisition pro- 
grams by entering the equation's coefficients. This could 
easily be implemented by prompting the operator for the 
equation's coefficients. The creation of a more automated 
method for base-line calibrations could provide a corrector 
curve and/or tape for direct application to the acquired 
navigation data. Automation of the current base-line cali- 
bration procedures and the ensuing computations that yield 
the calibration curve would enable the operator to sample 


more calibration data in the same amount of time. 
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3. Present Calibration Problems 
Some issues need to be addressed regarding present 
NOS methods of determining range error correctors for cor- 
responding signal strengths: 
The assumption of a best-fit linear curve being most 
representative of an MRS III's range error-signal 
strength relationship is not necessarily the most 


accurate one considering the variety of conditions 
and equipment configurations over which a system can 


Operate. Another choice is to use a separate cor- 
rector for each signal strength in the calibration 
curve. This approach assumes that a functional re- 


lationship exists between range error and signal 
strength and can be defined through curve-fitting 
procedures. This thesis discusses what degree 
polynomial best fits the January 1984 data after 
curve-fitting techniques were tried. 


The optimum base-line length for a base-line cali- 
bration has been revised several times by NOS. 

As new field phenomena are realized the recommended 
calibration length changes [Ref. 8]. When base 
lines as short as 100 meters are used, corrector 
anomalies occur at the higher signal strengths. 

Me converse situation Of observing range anomalies 
at low strengths for longer base lines limits the 
base-line length on the high end. The cause of 

the anomalies at close-in ranges is "the measured 
pulse amplitude riding on top of the previously 
decaying pulse." [Ref. 8] Anomalies occur at 
ranges less than 700 meters where long distances 
are being simulated with a variable attenuator. 

The dilemma is if too short a base-line length 

is selected, the resultant signal strength versus 
range corrector curve may not become exponential at 
low signal strengths and thus MASS cannot be determined. 
On the other hand, if too long a base-line length 
is selected, no correctors will be obtained at 
higher signal strengths and correctors for field 
data must be extrapolated from lower signal strength 
correctors [Ref. 9:p. 21. 


At present in NOS, MASS is a subjective estimate of 
"ο ροτς ος Inrlection (or signal strength) where 
the curve becomes exponential. A mathematical 

definition of MASS that can be used in future cali- 
bration software is desirable. Although the aim is 


ο]. 


to fit a continuous function to the calibration ons 
the curve is really not a discrete one in that non- 
integer signal strengths do not exist. In the case of 
the newer Falcon system, every signal strength may not 
be calibrated due to the expanded signal strength 
scale. The curve is thus empirically derived from 
those strengths that were calibrated. The MRS III 
signal strength scale is from 0 to about 25, while 

the Falcon scale range is from 0 to 100. 


Present base-line calibration procedures are not re- 
sponsive to observed anomalies in the calibration 
curve because a single mean corrector is selected 
for the entire curve. By using a functional form of 
the entire curve, anomalies could either be rejected 
through an outlier statistical test or become incor- 
porated into the correction process. The presence 
of these anomalies is confirmed by data obtained for 
this thesis. Anomalies are shown in the range-error 
versus signal-strength plots in Appendix B. 


This research should aid NOS in developing a calibration 
procedure that realistically predicts the range corrector 
for an observed range rate given the range and signal 


strength over a variety of system configurations. 
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IT. "MINI-RANGER III POSITIONING SYSTEM 


А SYSTEM COMPONENTS 
A good understanding of the components of the MRS III 


(Table III) and their operation will be helpful to the 


reader. 
TABLE III 
MINI-RANGER COMPONENTS AND SERIAL NUMBERS 

System Component Serial Number 
Range Console (RA-4), BO269 
RT ОПШЕ (КА-5), В1108 
Signal Strength Counter No. 1 WLO 

Signal Strength Counter No. 2 η 
Transponder Code 0 C1789 
Transponder Code l C1883 


l. ReccMver.Hwum5mittergeempby 


The term universal station refers to the assembly 
that houses either an RT unit or a reference station (Fig. 
2.1). A universal station contains a motherboard, power 
Supply, modulator, and the RF assembly. An RT unit is a 
universal station with a video interface assembly and a 
6-db antenna. The receiver-transmitter assembly (RT unit) 
consists of a radar transmitter, a radar receiver, anda 
video interface assembly. Due to the complexity of the 
electronics, technical descriptions of the system's compo- 
nents and their operation are taken directly from the MRS 
III Positioning System Maintenance Manual. 
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The encoded pulse modulation signal from the range con- 
sole passes through the video interface assembly to 
the transmitter. The resulting coded RF signal passes 
through the circulator to the antenna. At the same 
time, the transmitter causes a video signal to pass 
through the video interface assembly to the range con- 
sole where it becomes a start signal for the range 
counters. Meanwhile, the interrogation signal is re- 
Ceived, processed, and a reply is generated by the 
selected reference station assembly. The reply signal 
is received by the RT unit antenna and is passed 
through the circulator to the receiver. The resulting 
video then passes through the video interface assembly 
and the interconnect cable to the decoder assembly in 
the range console. There it is decoded and becomes 
the stop signal for the range counters. The receiver- 
transmitter operates from a +28 and -12 volt DC supply 
provided by the range console. The receive and trans- 
mit frequencies are factory preset to standard 
frequencies.  [Ref. 10:р. 7] 


2. Shore Reference Stations 

An RT unit is located on the survey launch and in- 
terrogates one or more reference stations located ashore. 

A reference station is a shore-based universal station that 
is located over a known geodetic control station. 

A reference station contains a l6-code decoder as- 
sembly and a 13-db antenna. The description of a reference 
station's operation is helpful in understanding possible 
error sources. An interrogation signal from the RT unit is 
coupled through the reference station receiver to the 16- 
code decoder assembly. If the proper first-to-second pulse 
and second-to-third pulse spacing is detected, the decoder 
will respond with a two-pulse signal that is coupled to the 
modulator, causing a reply to be transmitted back to the RT 
unit. The first-to-second pulse spacing is given in micro- 


seconds and corresponds to a given code value. The 
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second-to-third pulse spacing is either 50 or 54 microseconds. 
A reference station must be set to the code that has been se- 
lected on the range console. 

3. Antennae 

There are three different types of antennae for use 
with the MRS III (Figs. 2.2, 2.3, and 2.4). The antennae 
are classified by their purpose which is related to the out- 
put gain and radiation pattern characteristics. Elevation 
and azimuth patterns refer to the vertical and horizontal 
radiation angles of an antenna (Fig. 2.5). 

The medium-gain antenna (Fig. 2.2) has a gain of 13 
db and is the antenna most often used for hydrographic sur- 
vey operations and base-line calibrations. The horizontal 
azimuth pattern of the medium-gain antenna is 80 degrees and 
its> vertical elevation pattern is 15 degrees. 

The high-gain antenna (Fig. 2.3) has a gain of 19 db 

and offers the user a greater operating range. The high- 
gain antenna has a vertical pattern of only 5 degrees. 
The narrower vertical pattern focuses the radiated energy to 
increase the operating range of a shore station. The level- 
ing of a high-gain antenna is more critical then, due to its 
narrower vertical pattern. 

The omnidirectional antenna (Fig. 2.4) has a gain of 
6 db and is used on the survey vessel with the RT unit. It 
can also be used on shore reference stations located on is- 


lands or prominent headlands where 360-degree azimuth 
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coverage is desired. The vertical pattern is 15 degrees. 
An interesting fact about this type of antenna is its mis- 
leading name. "An omni antenna will receive signals from 
all directions. . . . The catch, however, is that it does 
not receive equally well from all directions. The omni an- 
tenna has a preferred direction." [Ref. 2:p. 8] 

Previous research has shown that the 90/270-degree 
directions are the preferred ones. Unfortunately, the  an- 
tenna is usually installed with the 0/180-degree axis point- 
ing toward the work area. (Ref. l10:p. 8] The January 1984 
calibrations were performed using the omnidirectional an- 
tenna with the 0/180-degree direction pointing toward the 
other antennae. 

4. Range Console 

The range console (Fig. 2.6) consists of a coder 
assembly, range control assembly, range counter display as- 
sembly, motherboard assembly, and a power supply assembly. 
The power supply assembly provides operating voltages for 
both the range console and the RT unit. A description of 
the operation of the range console follows. The coder as- 
sembly generates the transmitter pulse modulation signal. 
When a start signal is received from the receiver-transmitter 
assembly, the range counter starts to count. The reference 
station assembly reply video is decoded by the decoder and 
becomes the stop signal for the range counter. After five 


consecutive reference station assembly replies, the range 
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coter control commands the front panel range counter dis- 
play to refresh the range readout. The read information, in 
parallel binary-coded decimal (BCD) format, also is available 
at a rear panel connector for use by peripheral equipment 
such as printers or computers. The range console interro- 
gates each of the two reference station assemblies alternate- 
ly. Two thumbwheel switches on the front panel select two 
reference stations by code. Four distinct reference station 
assembly codes may be used, two at a time, in a standard MRS 
III. жапкорінбопте Со this alone selection from any of 16 codes, 
rather than four, to improve the versatility of the system 


Should other users be present in the work area. 


В. ЖЕС ІТ 5ҮСТЕМ ОРТТОМӘ 
1. Range Average 

An optional feature of the range console included in 
all NOS MRS III's is the range average switch. The range 
average switch performs internal averaging of the received 
range rates at one of five accuracy levels. This option 
helps the system produce more repeatable readings. 

A range average switch is located on the front panel 
of the console that is labeled in total turn-around cycles 
called PRF cycles. The five settings are l, 5, 20, 40, and 
75. "Іп the OFF position, the range console will take the 
normal 5 consecutive PRF cycles to produce one valid 
range reading. In the 75 position, the console takes 75 PRF 


cycles or 15 samples to automatically provide an average. 
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This is the same as the operator taking 15 readings from the 
front indicators and calculating the average. If the range 
console misses one reading in acquiring the average, the en- 
tire group of readings is nøt thrown out, only cHe Greup "or 
5 іп which the average was missed."  [Ref. 3:Section 5.2.2] 

Accuracy claims (by Motorola) for Έπε ΕΚ ΣΤΙΤ ασε 
based on a one-sigma standard deviation or probability of 
68.23 percent (Appendix A). The best results are obtained 
by using a single survey vessel and the 75 range average 
setting. "If more than one range console is used in the sys- 
tem, some erratic readings will result if the RANGE AVG 
switch is set above 20." [Ref. 3:Section 3.1.4] Higher 
settings require more update time. This is important if 
multiple users need to share shore stations. To determine 
whether range averaging is significant during a base-line 
calibration of a system, the January 1984 calibration data 
were acquired using settings of 20, 40, and 75. 

2. Signal Strength Indicator 

The signal strength indicator is mounted atop the 
range console (Fig. 2.7). The ranges shown on the console 
correspond to signal strength values displayed on the indi- 
cator. A brief explanation of the counter will aid the 
reader in understanding the range error/signal strength 
relationship. 

The video processor inside the RT unit converts the 


received pulse amplitude to a standard amplitude proportional 
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to the amplitude of the received pulse. This pulse is then 
used to start and stop a counter clock that generates a 
relative strength value in integer form [Ref. 3:Section 5.2.3]. 
Typically the signal strength will vary from 4 to 
25 for the MRS III. NOS counters are also equipped with audio 
alarms that are activated when signal strengths fall below a 
preset value. When the range console receives a signal that 
has undergone attenuation, it's degraded amplitude is con- 
verted to a relative value, displayed on the indicator. 
Signal strengths can be used with the alarm option 
to identify weak areas of coverage. NOS signal strengths 
are recorded directly on the data tape. The next step is to 
use these recorded signal strengths to extract a range-error 


corrector in nedr real time. 
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ІІІ. EXRERIMENT) DESIGNVAND DATA REDUCTION 


СО ЧОБ ТЕСТ Г У?ЕС 

The main objective of this thesis is to generate an MRS 
III Corrector curve that, when given a range rate and corre- 
sponding signal strength as input, will yield an accurate 
range corrector that can be applied to the range rate in 
near-real time. The corrector curve should be characteris- 
tic of all possible MRS III equipment configurations and con- 
sole settings. To obtain such a curve, the thesis addresses 
several related topics: 


Range Average was examined for its significance in de- 
fining the calibration curve. 


The significance of the three different types of antennae 
on the base-line calibration was investigated. 


The significance of base-line distance was investigated 
in hopes of finding an optimum base-line calibration 
length. 


Ameeaduipment configuration that is mest characteristic for 
an MRS III unit was identified. 
B. ASSUMPTIONS 
Several assumptions were made when performing the field 
work and analyzing the data: 
Range errors were not related to power supply voltages 
Since the units furnished by RAINIER were equipped 
with “cut-off" switches that deactivate the transponders 
when the voltage is less than 22 volts. Therefore, 
errors due to power fluctuations were assumed negligible. 
The horizontal and vertical pointings of the reference 


station antennae were considered to be negligible since 
BgenMeode/ reference stationawas pointed in the same 
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manner. The antennae were pointed visually for the three 
shorter base lines. For the two longer base lines, a 
radio was used to tell the reference station operator 
when optimum signal strengths were observed. 


The tellurometers that were used to measure the base 
lines were accurate to plus or minus 0.015 meters. 
Thus, the lines were measured to submeter accuracy 
so reference distances introduced no significant 
error into the data. 


Artificially attenuating the incoming signals to the 
RT unit adequately represented real conditions en- 
countered in the field. Though this assumption is 

not necessarily valid, the artificial attenuator 
offers the only practical means of simulating the 
attenuation caused by signals traveling long distances. 


Internal drift of the MRSIII's duringwbthescalubnhtienus 
was assumed to be negligible. Each day a warm-up time 
of 30 minutes was observed before the first calibration. 
Throughout the day, the units were not powered down and 
antenna changes were affected by momentarily halting 

the interrogation process with the “Hold Display" button 
on the front of the range console. Morning warm-up 
times were within the prescribed time as recommended in 
the NOS Operation and Installation Manual. Thus, range 
errors due to temperature changes within the magnetron 
of the MRS III system were also considered insignifi- 
cant since ambient temperatures during the January 
calibration were invariably within two to three degrees. 
(Ref. 3] 


C. DATA ACQUISITION AND PROCESSING 

The calibration field work for the thesis was performed 
January 27 to 30, 1984, across Puget Sound and Shilshole Bay, 
Washington (Fig. 3.1). Base-line calibrations were performed 
over five different distances uSing the station at West Point 
Lighthouse as a pivot point (Table IV). 

l. Tellurometer Measurements 


Prior to beginning the actual static calibrations, 


each base line was measured to submeter accuracy uSing a 
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Figure 3.1 Project Area and Station Locations 
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TABLE IV 


BASE-LINE ENDPOINTS AND RESPECTIVE LENGTHS 


(From West Point Lighthouse) 


Base Line End Point Length(m) 
A Sewage Gate 1061115 
B Shilshole Marina 2417.48 
e Meadow Point 4083.10 
D Point Monroe 7016. Т 
P: Point Jefferson 9861.08 


Tellurometer model CA-1000 electronic distance measuring 


instrument. 


Since the MRS III range rates were slope distances, 
no slope to horizontal distance computations werè necessary. 
However, due to the tellurometers having a preset value for 
the refractive index in air of radio waves equal to 1.000325, 
weather data were taken concurrently to correct this fixed 
value for variations in wet- and dry-bulb temperatures and 
for barometric pressure. These same weather data were input 
into a Hewlett-Packard 9815S calculator. Meteorological 


correctors were performed using an NOS geodetic program 


(EDM-03/821201, option A). 


The Tellurometer observation forms (NOAA Form 76- 


61A) for the five base lines include the output tape gener- 


ated by the HP 9815S calculator (Appendix C ). 


These submeter results were used as reference dis- 


tances for the five base lines and subsequently in the de- 


termination of average range errors to 0.1 meter. 
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2. Mini-Ranger Data Acquisition 


Base-line calibration procedures are documented in 
the PMC OPORDER [Ref. 5:Sections M-3=through Μ-5]. The 
range console and two signal strength counters were placed 
in the rear opening of the rental hatchback to shield the 
recorder from strong winds. Two reference stations were set 
to codes l and 0, respectively, and were positioned at the 
designated endpoint (Table IV) where an operator aligned and 
pointed the reference stations at the station near West Point 
Lighthouse and changed antennae when requested. 

For each signal strength setting, ten range rates 
were recorded. The RT unit was artifically attenuated and 
the resulting range recorded for each corresponding signal 
strength. After calibrating a code for three range average 
settings--20, 40, and 75--the medium-gain antenna was re- 
moved and replaced with a high-gain or omnidirectional 
antenna. 

All three antenna types were calibrated for each of 
two codes at three range average settings for each of five 
base lines. The total number of calibration "sets" was 3 x 
5 х 3 x 2 = 90. The beginning time of each individual cali- 
bration was recorded and was included in the header record 
for each data file. 

The range console displays A and B show independent 
rates. This allows the observer to sample two independent 
rates at once. For each calibration set, ten readings were 


obtained by pressing the "Hold Display" button five times. 
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The ten readings were used to compute the mean range error 
for each signal strength. 

3.  Mini-Ranger Data Processing 

| Data from the January 1984 calibrations were manu- 
ally logged into data files on the Naval Postgraduate School's 
IBM 3033 computer. Ten readings per signal strength were re- 
corded for each code (codes 0 and 1). Each data file con- 
tained a header record that identified the base line, code, 
range average setting, antenna type, and time of calibration. 

The data were printed ina readable form and checked 
for correctness. The data were also corrected for a system- 
atic offset of 16 meters that was not adjusted or zeroed. 
Instead of adjusting the console, the l6-meter offset was 
applied to all the data. 

After verifying the correctness of the data, a 
FORTRAN/DISSPLA [Ref. 11] computer program STAT (Appendix D) 
computed mean rates, mean range errors, and variances for 
both codes. 

A "calibration set" is herein defined as a distinct 
set of signal strengths and corresponding range errors for a 
specific configuration of code (0 or 1), antenna type, range 
average setting, and base line. 

Program SSPLOT (Appendix D) plotted mean range error 
versus signal strength and standard deviation for each mean 
error (Appendix B). Standard deviation is displayed as a 


vertical bar for its respective range error. 
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Plots of range error versus signal strength were 
generated to assess the ОТО ОЕ Салса of range average and 
m enna type to the- -overall calibration curve for each code. 
Program SCATTER (Appendix D) generated scatter diagrams 
(Appendix B) of mean range errors for each of the two codes. 

Program WALK (Appendix D) generated plots of mean 
range error versus base-line distance using signal strengths 
as markers. These plots (Figs. 3.2 and ο ο) Show that elower 
signal strengths consistently separate themselves from the 
others. | 

E unctional Curve Fitting 

Curve fitting is a technique that employs a least 
squares criterion to decere unique regression curve of 
any degree up to a n minus l, where n is the number of paired 
data points. The fitted Curve rs a functional form or equa- | 
tion from which range correctors can be predicted in real 
time. SDS III would be capable of referencing electronic 
Correctors by two methods. One method uses a table of 
stored data points in an array; it is the simplest but would 
require some type of interpolation for intermediate signal 
strengths. A second method uses the functional form of the 
calibration curve to extract range errors for a given signal 
Strength. 

The main advantages of the functional approach are 


Τετ computer computations and. the accounting 
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for intermediate signal strengths in a nonlinear, more accu- 
rate manner. The major reason to limit the time spent cali- 
brating electronic control systems is to maximize data 
acquisition time. NOS's desire is to minimize the calibra- 
tion time by using faster, noncritical daily system checks. 
If the calibration procedures were automated, a greater num- 
ber of data could be sampled and several base lines or equip- 
ment configurations could be investigated. The time saved 

by letting a mini-computer compute the statistics and re- 
gression curves directly from a floppy disk containing the 
raw calibration data would enable the operator to obtain more 
calibration data in much less time. A curve-fitting program 
POLY (Appendix D), obtained from Gerald's book on numerical 
analysis, was applied to the data [Ref. 12:p. 492]. 

The Gerald program computes coefficients for up to a 
seventh-degree polynomial regression curve, though only fifth- 
degree and lower were attempted on the January calibration 
data. Both the beta function (variance) and the graphical 
presentation suggest that a third-degree fit was the highest 
polynomial fit obtainable. Graphs were made that show re- 
gression curves for powers of one through three along with 
the actual mean range errors that were used to produce the 
coefficients defining the regression curves. These graphs 
are discussed in Chapter 4. Since program POLY does not 
print the regression coefficients, a Hewlett-Packard 41C 


statistical program (SIGMAPOLYC) was used to compute the 
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coefficients for a cubic regression curve for each base line 
for the two codes. А11 data and programs used in the thesis 


have been recorded on 9-track magnetic tape. 
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IV. OBSERVATIONS AND CONCLUSIONS 


A. OBSERVATIONS 
1. Scatter Plots 

Scatter plots (Appendix B) show range ἐν τῶν. versus 
signal strength for each calibration set. By inspection, 
base lines D and E contributed very little to the scatter 
plot for both codes. Range error differences of up to 9 
meters were observed at the two lowest signal strengths. 
Base line B exhibited the least sample deviation and covered 
the greatest range of signal strengths. Though fewer data 
points were available for the two longer base lines D and E, 
the general shapes of these curves were similar to the three 
shorter base lines. 

Data anomalies occurred for base line D code 0 (Ap- 
pendix B) between signal strengths five and eight. A maxi- 
mum difference of 5 meters was observed between the 
high-gain antenna curves and the other two antenna types. 
The agreement between the three antenna types was typically 
1 to 3 meters and is best seen on the mean range error 
plots. 

2. Mean Range Error Plots 

Mean range error plots showed differences in mean 

range error due to the range average setting were not signif- 


icant across base lines or antenna types. A visual 
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inspection of the data showed excellent agreement between 
different range average settings. 

High-gain antenna data showed the least sample de- 
viation in a consistent ov This fact and the great 
range of signal strengths over which one can calibrate sug- 
gest that calibrating with the high-gain antenna produces a 
calibration curve representative of all three antenna types. 
An anomaly in the high-gain antenna data occurred in base 
line A (around 1,000 meters) for code 0 (Appendix B). The 
plot shows an upward shift of the curve to range errors of 
3 to 5 meters higher than the data from the other 
antennae. The same behavior was not observed with code l. 
Inspection of the mean range error versus base-line distance 
Blots (Fags. 3.2 and 3.3) point tova MASS of signal strength 
five for all the plots. 

3. Range Error Versus Base-Line Distance Plots 

The range error versus base-line distance plots 
(Figs. 3.2 and 3.3) show a dip where range error decreased 
when signal strength increased from four to six and then in- 
creased for signal strengths seven and above. This charac- 
teristic dip near signal strengths four and five was noted 
for all equipment configurations and distances. Signal 
strength four was separated from the remaining signal 
strengths by a constant positive value of 3 to 4 
meters for each base line. The other range errors (with 


signal strength markers) were grouped close together. The 
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ο RE dip mentioned above was not as pronounced on 
the fitted regression curves--especially the cubic fit. 

As anomaly occurs on the base line B range error 
versus base-line distance plot. Base line B was partially 
obstructed by a pier 200 yards forward of the West Point 
endpoint. The above mentioned range errors may have been 
the result of additional attenuation from the pier. 

4.  Curve-Fitting Results 

Regression coefficients were computed for degrees 
one through four (Table V). Variance was also computed in 
the program for each degree of fitted polynomial and was 
used as a numerical indicator of a curve's "goodness of fit" 


(Table VI). The regression curves were plotted (Figs. 4.1 and 


TABLE V 


REGRESSTONICOEFF TCI ENIS 


(y = a + bx + E” + а + EE 
Degree of Coefficients 
Polynomial Variance a b С d e 

Code 0: 

| -0.88 0.88 0.00 ο ου 0-08 

2 4.54 PPS ο ος ο ο μυ WU 

3 INA. 38 -5.66 0.59 =0%02 ο 

4 1754 -5.66 0.535 -0302 0.00 
Code 1: 

1 On OT 205192 0.00 0 001 ОТОО 

2 ΕΘΞ -0.70 ors 52. 202200 

3 1.324 -4.06 0538 00 О 

4 1354 -4.06 0.39 2001 000 
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TABLE VI 


COPPETICTENTSIFOR CUBIC REGRESSION CURVES 


(y =a + bx + eeu + dx?) 


Base Line  VARIANCE a b © а 
Code 0: 
А ο. 53 18.56 -7.47 0.84 -0.03 
B 0.64 18.14 -5.36 ο. οι -0.02 
С 0.69 29 02 -7.86 032 -0.03 
D 9755 53.44 -21.08 2.60 -0.10 
Е 2 247559 150% 2.27 ος 
Сойе 1: 
А (ess 9.11 -2. 94 0.28 -0.01 
B 0:52 M5 οἱ -4.22 0.40 -0.01 
С ОС 25,49 -8.97 0.95 -0.03 
D ο οὐ 66. 84 -27.09 5.5 =0 115 
E 0.80 42.67 = 77 2.63 -0.12 


4.2) with the mean range error/signal strength pairs used as 
input for creating each higher degree fit. Gerald's book 
[Ref. 12:p. 437], describes the variance as an indicator for 
deciding what degree of polynomial should be used. As higher 
degree polynomials are used, the deviations of the points 
from the fitted curve will decrease until, when the degree 

of the polynomial n equals n minus 1, there is an exact 
match (assuming no duplicate data at the same x-value). 
Variances for the four higher polynomial regression curves 


me ated tnat the cubic fit was the best. Опе 
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increases the degree of approximating polynomial so long as 
there is a decrease in the variance. Both codes had lowest 
variances for the third-degree curve. Knowing this, cubic 
regression coefficients were computed for each base line for 
the two codes (Table VII). 

The coefficients for the base lines were computed us- 
ing a Hewlett-Packard 41C calculator and regression program 
( SIGMAPOLYC). 

The good agreement between cubic fits at base line B 
for each code and their respective overall curves suggests 
that a distance of about 2,500 meters is an optimum range over 
which to calibrate (using endpoints at sea level). The cubic 
regression curves for this distance adequately depict range 
error versus signal strength for all five base-line distances 


and three types of antennae. 


B. CONCLUSIONS 

The optimum base-line length for obtaining a representative 
calibration with both calibration endpoints at sea level was 
Closest to base line B (2,500 meters) for both codes 0 and 1. 
At present, NOS requires the base-line length to be between 
700 and 3,000 meters. 

The high-gain type antenna varied slightly more than the 
other two antenna types.  High-gain calibration errors agreed 
to within l to 2 meters with the medium-gain and omnidirec- 


tional antenna results in all but 6 of the 90 calibration sets. 
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ШВЕИ а 


COMPUTED VARIANCES FOR HIGHER DEGREE CURVES 


Polynomial Degree Variance 
Code 0: 
Ji IG 
2 05 
3 0.43 
4 0.43 
Code 1: 
ii EO 
多 0.. 87 
3 0.54 
4 0.54 


The fitted polynomials (Figs. 4.1 and 4.2) had maximum 
discrepancies in range error of 1.5 and 2.0 meters for codes 
0 and 1, respectively. To improve the quality of fit, one 
could increase the number of data points. 

Another possibility is to resolve a calibration curve into 
exponential and linear components and fitting each component 
separately. The expanded scale of the Falcon 484 offers a 
Greater range of signal strengths over which to calibrate. 

The improved mean values would improve the fit of the function 
to the observed data. 

It may be prudent to set operational distance limits for 
the medium-gain and omnidirectional antennae. Both antennae 
were lacking in data points at the higher signal strengths 
for calibrations at distances greater than base line B (2,500 


meters). Likewise, a minimum distance may be needed for 
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high-gain antennae as shown by the code 0, base line A range- 
error plots (Appendix B). The anomalies suggest that cali- 
brating with a high-gain antenna at too close a distance 
(perhaps less than 1,000 meters) may result in erroneous 


correctors. 


C. | RECOMMENDATIONS 

Study the concept of rejecting data below MASS and fitting 
the remaining data to a linear curve. | 

Investigate spline curve fitting or the resolving of a 
calibration curve into linear and exponential components be- 
fore deciding upon future calibration procedures. 

Give consideration to mathematically determining the MASS 
of a curve prior to fitting. Further statistical tests 
could then be applied to data below the MASS to determine 
whether it should be rejected. 

Plot mean range error versus signal strength for each 
base line and for the cumulative mean range errors for each 
of the two codes, thus showing which base line's data in- 
fluences the overall mean curve the most. A comparison 
should be made using a single, mean corrector MRS III code, 
as is presently done, with the predictive results obtained 
from curve fitting. The comparison should consider the 
allowable positioning errors for different survey scales. 

Review present base-line calibration procedures with 
respect to the findings of this thesis. Specifically, ral: 
uate the optimum distance of 2,500 meters with a range 


average setting of 75 and a high-gain antenna. 
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APPENDIX À 
MINI-RANGER III SYSTEM SPECIFICATIONS 


Overall System Power Requirements 
Range Console 
Input Voltage 105 to 125 volts AC 
12 το πε 
Power Consumption ` 60 watts AC 
40 watts DC 


RE Upit 





Input voltage 22 to 32 volts DS 
Power Consumption (at 27 VDC) 
No interrogation (standby) 9 watts 
Normal interrogation l6 watts 
Transponder 
ο. ου voles DC 
Power Consumption (at 27 VDC) 
No Interrogation (standby) 8 watts 


Normal Interrogation 15 watts 
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lransmission Frequencies 


P 
G 
5 
H 
ct 





Transmit 5570 MEZ 


Receive 5480 MHz 


Shore Iransponder 





Transmit 5480 MHz 


Receive 5570SMAZ 


Number of Users 


With range averaging and multi-user options installed, 
Motorola recommends no more than five range consoles operate 
in a project area. Although the consoles may have different 
codes, they will still lock each other out when transmitting 


to a shore station. 


Accuracy of the System 


The following figures are based on one-sigma standard 
deviation. This means that 68.23 percent of all readings 


will fall within the quoted accuracy. 


Range Average Setting Accuracy (in meters) 
й 3. 00 
5 212 
20 ТЕБ 
40 1. 06 
75 0.78 
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Antennae and Beam Characteristics 


Receiyer-Transmitter (RT unit) 


Gain 6 db 
Horizontal 360 degrees (omnidirectional) 
Vertical 25 degrees 


Shore Transponder 





Medium-Gain Antenna 
Gain оосар 
Horizontal 80 degrees 
Vertical 15 degrees 
High-Gain Antenna 
Gain натар 
Horizontal 80 degrees 
Vertical S degrees 
Omnidirectional Antenna 


Same characteristics as the RT unit antenna. 


Range of System Operation 


The system is line-of-sight. Though Motorola's opera- 
tion manual claims a maximum range of 37,000 meters, this 
distance is undoubtedly achieved by increasing the elevation 


of either or both the RT unit and shore transponder 


antennae. 
Antenna Type Maximum Range (in meters) 
Medium-Gain 227149 
High-Gain 64,360 
Omnidirectional unknown but less than Medium-Gain. 
Motorola Options used by NOS 
Option Description Option Number 
Signal Strength 2 
Range Averaging 4 
Multiuser 5 
Byte Test Module 11 
30-Minute Standby 15 
16-Code Option 46, 47 
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APPENDIX B 


RANGE ERROR AND STATISTICAL PLOTS 


The following graphs are (1) range error versus signal 
strength (including standard deviations) for each calibration 
set (pages 72 through 101) and (2) scatter diagrams of mean 
mg qe rror “(for Gach signal strength) versus signal strength 


for each of two MRS III codes (pages 102 through 113). 
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APPENDIX C 


TELLUROMETER DATA 


The following observation forms show measurements of the 
five base lines obtained with the Tellurometers. All dis- 
tances were corrected for meteorological effects as shown on 


the attached HP 9815S calculator tapes. 
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NOAA =ORM 76—-61A 
(12—73) 


U.S. DEPARTMENT OF COMMERCE 
NATIONAL OCEANIC AND ATMOSPHERIC ADMINISTRATION 


TELLUROMETER OBSERVATIONS 
MODEL CA1000 


μα. Γ᾿ Ῥῴξετοονο |29 Jan 64 
STATION PHEASANT INST, NUMBER HEIGHT OF INST. 
[Bos ΚΒ 


WEST POINT E 


PROJECT NUMBER CHIEF OF PARTY OBSERVER . 


METEOROLOGICAL DATA FIELD READINGS 


ool ano Be 
Језа 
БЕСТЕН 


Siope distance 











т 
«4 
Fe ace эса Ал — — tsa "t t 
PREESS= zbh,.208lH, 
| NM - nm i MEAN SETS 1&2 — > 4 | 117 7Ет 
ШЕШЕ = 4s, ο δα г 
| li Fi E= w. ΒΡΕ 
do 
i ПЕЦ 
И τή -тыт- o nau hi 
REM. - „Гл Hat 2 aM METERS FEET X 0.3048006096 
Peele Sie Bee ы 
рен ieee = = mia 
| -— τε = t = - - - 132 | 
DISTANCE (Meters) | | a | ° ZAG үя 
τας Ὁ e - —— M ———— ων ο ' τ 
Тез: t — EE T τα 
| [PO ss poe ed l | 
E iR αἰ ι METEOROLOGICAL | 2. 
C.F H = - nes ғ. ' m 
| z Μος pue ar = 31 CORRECTION | С) оо 
- 1 ањ i ос. ЕКЕШ Si в = 
= -3.B8c5H | 4) 
NN cu 2 FIXEO INSTRUMENT 
σε РИТ - “IF = | OOES | 
- - κ: ΠΝ CORRECTION ° | 
ЗЕ ш ц, : | =; = 
CORRECTED | ( / 76 
SLOPE DISTANCE | “ο (o e 


ΠΡ 


NOAA FORM 76-61A U.S. DEPARTMENT OF COMMERCE 


(12-73) NATIONAL OCEANIC ANO ATMOSPHERIC AOMINISTRATION 






TELLUROMETER OBSERVATIONS 
MODEL CA1000 


RE s 
STATIO SHILSHOLE PER INST. NUMBER HEIGHT OF INST, 
Eese КШ 
Шат |С» _ 
— сты: 


METEOROLOGICAL DATA FIELD READINGS 


| SET | SET-2 
T T 


јео зря 
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Program 


METOBS 


RAWRITE 


STA; 


SSPLOT 


WALK 


P@l YF 1 T 


POL v 


BASELINE 


APPENDIX D 


COMPUTER PROGRAMS 


Description 


Uses DISSPLA subroutines to generate plots 
of temperature and barometric pressure 
versus time. 


Reads Mini-Ranger rates from a file; 
corrects the rates for a 16-meter offset: 
and writes the corrected rates, header 
information, and signal strengths to 
another file. 


Lists the means and variances of both 
observed rates and computed range errors 
of Mini-Ranger calibration data. 


Uses DISSPLA subroutines to generate plots 
of mean range error (Cand standard 
deviation) versus signal strength. 


Uses DISSPLA subroutines to generate plots 
of mean range error versus base-line 
distance. 


Fits polynomial curves for up to degree 
seven to the mean range error values. 


Uses DISSPLA subroutines to generate plots 
of range errors predicted from polynomial 
regression curves. 


Uses DISSPLA subroutines to generate plots 


of range error versus signal strength for 
each base line. 
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X e e ek kh he oc Ck kk kh o C e GG 


PROGRAM: METEOROLOGICAL OBSERVATIONS VS TIME (IN HRS.) 

AUIHOR: ЕКШСЕ ΕΤ ΠΙΤΕΡ ΔΕΟ 

DATE : APRIL 1984 

PROGRAM DESCRIPTION: THIS FORTRAN PROGRAM USES DISSPLA 
SUDEDUMNIESSESFEOTGENERBTE PLOTS. OF 
TENEERATURE AlID BAROMEGRIG PRESSURE 
VERSUS HOURS: 

ДИРОТ: INPUT COMES FROM DATA SIATENENTS. 


OUTPUT: ΤΗΣ ΕΙΡ Ισ COT CANSEE “DIRECTED IO ONE OF THREE 
БЕШШІЕРІ: VERomibewocc2A, VEKSATEC O900A, OR A 
ТЕКТКОМІХ 4631 HARDWARE DEVICE. 


Ἂ k & k k k k k k k ck ok X; k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k kk kk k k k k k k k k k k k k k kk k k k 


ΩΩωΩωωώωώωώώώώώώώώωώώώ 


REAL TEMP,PRESS,TINE 
DIMENSION TE!NP(50),PRESS(50), IPAK(100) , TINE(50) 


X k k k k < + k < K Jan. 27 DATA 

DATA О. 20700: L000 Sirs 00), foe e 72900... 2200. / 

DATA TEMP/47.0 ‚48. 0, 49. 0550. EN 51.0 55). 0, ο @ 

DAIA PRESS/30. 43, 30.40, 30. 42, 30. an ‚30. 40, 30.44,30.44/ 
k**X*kk**kk JAN.28 DATA 

ο ο τής οσους. 10ου. 1300. , 1600. 1909. 2200. / 

БАТА ТВЫР/49.0» 4870, 5250 7965 0250.0790-. 0 

ГАТА РКЕ55/30.48,30.46,30.44,30.38,30.40,30.40/ 
kXkkkkkkkk JAN. 29 DATA 

DATA TIINE/0400770700 999000. ου ^ 10007 72200* / 

DATA TEUP/47.0, 46. ο. 41. 0,47. D 46. ο 44.0 

DATA PRESS/30. 42, 30.42, 30. 42, 30. FOr SOS 5e30235/ 
AKAKKKARKAK JAN.30 DATA 

wae T MVE /O400"% 1000. 1308, 16059. ΒΟΟ,,.2200. / 

DATA TEIP/43.0,44.0,44.0,43.0,44.0,42.0/ 

eee ess /30.32,30.30,30.26,30.22 720.18 ,30.17/ 
kXkkkkkkkk INITIALIZE PLOTTER AND SET PAGE SIZE. 

САРЫ РЕШЕШ 72 6) 

CALL СОПРАЗ 

CALL πο το 0 ,0,0) 

CALL TEROIS 

CALL PAGE (11.,8.5) 
大 大 大 大 大 大 大 大 大 大 ]1AKE SS AXIS IN INTEGERS. 

CALL AINTAX 


AAKKAKAKKK ORIENT Y AXIS NUMBERS HORIZONTALLY. 


nananana 


ANA AAA AA AAAAAAAAAAAAAAAAAAAAANAQA 
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NANA aaqa (CG) G 使 MAM асай AAA 


ana 


Оз Оо СО а) 


өзе 


CALL YAXANG (0.0) 
kX******** SELECT TOTAL PLOT AREA. 
CALL AREA2D(8.,5.) 
kkkkkkkkkk SCALE AND LABEL AXES WITH SELF-COUNTING OPTION. 
CALL YNAME('TEMPERATURE (DEGREES F)$',100) 
CALL XNAME m 27, HOURSS',100) 
CALL BASALF('STAND' ) 
kkkkkkkkkk CREATE THE HEADING AND SET THE PRINT TYPE. 
CALL SWISSM 
CALL HEADIN ('METEOROLOGICAL DATA AT WEST POINTS',100,1.5,1) 
call duplx 
AkKAKKKKAKK DEFINE BOTH AXES AND FRAME THE SUBPLOT AREA. 
CALL GRAF(0200. ,200.0,2400. ,40.0,2.0,54.0) 
CALL ТНКЕВМ(.02) 
CALL FRAME 
大 大 大 大 大 大 大 大 大 大 SET UP LEGEND DATA AND TEXT FOR LEGEND NAME. 
CALL LINESP (2.0) 
CALL LINES рсе: ые. ш 
CALL LINES ('PRESSURES',IPAK,2) 
大 大 大 大 大 火 大 大 大大 FIND SIZE FOR LEGEND. 
ИЙ ОИ m 
YW=YLEGND(IPAK, 2 
CALL MYLEGN('OBSERVATION TYPE:S' ,100) 
AkKAKKKAKK DEFINE CURVE TYPES AND EXECUTE THE CURVE PLOT. 
CALL LEGLIN 
CALL DOT 
CALL CURVE(TIME,TEMP,7,1) 
大 大 大 大 大 大 大 大 大 大 SET UP SECONDARY Y-AXIS AT RIGHT END. 
CALL YGRAXS(30.00,0.05,30.60,5.,' PRESSURE (IN.)$',-100,8.0,0.0) 
CALL DASH 
CALL CURVE(TIME,PRESS,7,1) 
AX**XxX****k DRAW THE LEGEND. 
CALL LEGEND(IPAK,2,4.50,1.00) 
大 大 大 大 大 大大 大 大 大 CLOSE OUT METAFILE AND RETURN TO LEVEL ZERO. 
CALL ENDPL(0) 
CALL DONEPL 


STOF 
END 


eZ 


AAAAANAAAANANANANIN 


ОО О О 


AAAAAAAAANAAN 


大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 


PROGRAH: RAWRITE 

BESCR EE TLION: THIS РВОСВАН”ВЕАрС MINT-RANGER RATES FROH À FILE 
BUDAGORNECRSSTHENMFORENEIO HETER OFFSET. THE 
IOO JITHEN WBERBESOSTHESCORRECTED RATES, HEADER 
INFORMATION, AND SIGNAL STRENGTHS TO ANOTHER FILE. 

DATE. HAY 20. 1984 

AUTHOR: BRUGE F. HNPUARD 

INPUT/OUTPUT: BOTH USE À FILEDEF FOR READING AND WRITING. 

SUBROUTINES: NONE 


大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 


CHARACTER BL,A,B,C,D,E 

DATA A,B,C,D,E/'À','B' ,'C' ,"D! ,'E'/ 

INTEGER I,J,N,SS,RÁT,RATE,TIHME,RAVE,CODE,DATE,ANT,BASEL 

DIMENSION BÁSEL(1500] ANT (1500). CODE (1500), TIHE(1500),RAVE(1500), 
1 SS(1500),RAT(1500,10 7, RATE(1500,10), 

1 LENGTH(1500),DATE(1500) 

REAL LENGTH 


大 大 大 大 类 大 大 六 六 大 六 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大大 大 大 大 大 大 大 


PEGIN READING THE DATA FROM À FILEDEE INTO ARRAYS. 


类 类 类 六 六 六 六 六 大 类 大 文大 六 类 类 大 类 六 大 大 大 类 类 大 大 类 类 大 六 太太 大 大 大 六 大 大 大 类 类 大 大 大 大 大 大 类 类 大 六 六 大 大 大 大 类 类 大 大 大 类 大 大 大 大 
I=0 
7 I-I-*l 
READ(4,201) SS(I),(RAT(I, | 
1 RAVE(I), TIME(I 
Do TON"J-l1,10 


RATE(I,J)sRAT(I,J)-16 
100 CONTINUE 
IF(SS(I).EQ.99) GO TO 50 
Go To 7 
50 CONTINUE 


N=I 
大 大 大 六 大 大 大 六 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 六 大 大 大 


AT THIS POINT, THE DATA HAS ALL BEEN STORED IN THE ARRAYS SS(I), 
RATE(I,J),BASÉL(I),CODE(I),ANT(I),RAVE(I), AND TIME(I). 


六 大 大 大 大 六 大 大 类 大 类 大 大 大 六 大 大 大 大 大大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


太 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 六 大 大 大 六 


''' ΤΕ ΓΕΝΘΤΗ,ΒΙ, AND DATE AND BEGIN NO ОШТ ТНЕ DATA. 


大 六 大 大 大 大 大 大 大 类 大 类 大大 大 六 大 大大 大 类 大 大 类 大大 大 类 大 大 大 大 类 大 大 大 类 大 大 大 大 大 大 大 大 大 类 大 大 大 大 大 类 大 大 大 大 大 大 大 大 六 大 大 大 大 


Do 760 1-1,N 
IF(SS(I).EQ. 99) GO TO 500 
IF(BASEL(I). Go ТО 5] 
IF(BASEL СО ЛО 52 
IF(BASEL GO πο 53 
IF(BASEL(I). GG NO 54 
IF(BASEL(I).EO.5) GO TO 55 
51 LENGTH(I)21061.24 
BL-A 
ПЕШЕ(ТІ)-30 
GO TO 59 
52 LENGTH(I)-22417.48 


L-B 
DATE(I)-29 
ας. τὸ ο” 
93 LENGTH(I)=4083.10 
BL=C 
DATE(I)=27 


GO TO 59 
54 ІЕМСТН(І)-7016.77 


I-H P FA FA HH 
[rj 
tO O 
fS CN I 


T23 


ὦ ὦ ὦ С посао о ANA 


GQ O GO OOO 


BL=D 
DATE (T)= 28 
Comme 

55 PENCTH(I)= 9861.08 
BL-E 


DATE (I )=28 


CO IO 59 
大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 炎 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


CHECK FOR A NEW CALIBRATION SET (HEADER 
大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 类 大 大 
59 RUE .EQ.1) GO TO 65 

60 IF ο .LE.SS(I-1)) GO TO 80 

65 WRITE(6,202)BL,LENGTH(I),DATE(I),TIME(I),CODE(I),RAVE(I) 


大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 火 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


大 大 大 大 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


TEST FOR THE ANTENNA TYPE TO BE USED IN THE HEADER. 

大 大 大 类 大 大 类 大 大 类 大 大 类 大 大 大 大 大 大 大 类 大 类 大 大 大 类 大 大 类 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 六 大 类 大 大 大 大 大 大 大 大 大 大 大 大 
IF(ANT(I).EQ.1 
IF(ANT(I WRITE(6,204 


MITE 8 204) 
WRITES, 205 


.EQ.2 
IF(ANT(I 


οσο 3 
WRITE (6: ‚206 


WRITE 


i) 
大 大 次 大 大 大 大 ТАП 大 类 类 大 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


BEGIN WRITING OUT THE DATA ITSELF. 
大 大 大 大 大 六 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 


80 CONTINUE 
WRITE(6,208) SS(I),(RATE(I,J),J21,10) 


160 CONTINUE 


kc heck e ehe hh hee he he hh e ehe ehe he Ahhh ehh hh Rh hk 
FORMAT STATEMENTS: 
大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 
201 τ. 2Х,10(14,1Х), 11, ρω ЕЗ” 
202 που / 2x, ' BASE LINEA, 2X, ' PENGON = „,Е7.2,1Х, METERS’ , 322 
1 ΠΠ 27 А, 19841, 23 14 ÍX,'HOURS',/,18X,'CODE: ',Il,4X 


1 'RANGE AVERAGE SETTING: p 
203 FORAT, ' ANTENNA TYPE: MEDIUM GAIN!) 
204 FORHAT(22X,'ANTENNA ТҮРЕ: OHNIDIRECTIONAL!) 


205“FOBHAT ЛЛК 


206 FORHAT /,3X,'SIGNAL' 


207 HERHATUZXS 


'ANTENNA TYPE: HIGH GAIN!) 
,15X,'OBSERVED MINI-RANGER RATES!) 


ο ΕΕ σα ο... X,'2!.,4A,'3 ал "Т 


1 45,17! 45,181 45,191 45 (ШӘ 
208 FORMAT(4X,12,5X,10(14,1X)) 
500 CONTINUE 


ГЕО. 


STOP 
END 
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οκΧΧΑΧΧΧΧΑΧΧΑΧΛΑΧΧΧΑΧΧΑΧΧΧΧΧΧΧΧΑΧΧΧΧΑΧΧΑΧΑΧΧΧΑΧΧΑΧΑΧΑΧΑΧΑΧΑΧΧΑΧΧΧΧΧΧΧΧΧΧΧΧΧΑΧΧΧΧΑΧἍΧἍἌἍἍ 


AAAAAAAAAANAQAQ 


QAAAAQ 


AAAAAAAAAANQA 


PROGRAM: SIAT 

DESCRIPTION: THIS PROGRAM LISTS THE MEANS AND VARIANCES OF BOTH 
OBSERVED RATES AND COMPUTED RANGE ERRORS OF THESIS 
MINI-RANGER CALIBRATION DATA. 

DATE: MAY 20, 1984 

AUTHOR: BRUCE F. HILLARD 

INPUT/OUTPUT: BOTH USE A FILEDEF FOR READING AND WRITING. 

SUBROUTINES: NONE 


大 大 大 大 大大 大 大 大 大 大 炎炎 大 炎炎 大 炎炎 大 大 炎 大 大 火炎 大 大 大大 六 大 大 大 大 大 大 大 大 大 大 交大 大 大大 大 类 大 大 大 大大 大 炎炎 大 炎 大 大大 大大 大 大 大 


CHARACTER BL,A,B,C,D,E 

Dama A,B,C,D,E/°A', "Bo ЕТУ 

INTEGER I,J,N,SS,RAT,RATE, TIME, RAVE, CODE, DATE, ANT, BASEL 
DIMENSION BASEL (1500) ANT(1500$ , CODE(1500) , TIME(1500) ,RAVE(1500), 
1 $S(1500) ,RAT(1500,10) ,RATE(1500, 10), ҚЕМІЕ |1300). ,SUMM(1500), 
1 LENGTH(1500), DATE 1500), /XRERR(1500) , VRATE(1500 

1 VRERR(1500),SUMMM(1500),RERR(1500,10) ,RATSUM(1500) , SUM(1500) 
REAL LENGTH, RERR,RATSUM,XRATE,SUM,SUMM, SUMMM , XRERR, VRATE , VRERR 


ἍχχΧχχχΧΧΧΧΧΧΧΧΧΑΧΧΧΧΧΑΧΧΧΧΧΧΧΑΧΧΑΧΧΑΧΧΧΧΧΧΧΧΧΑΧΧΧΧΧΧΑΧΧΧΧΧΑΧΧΧΑΧΧΧΑΧΧΧΧΧΧΧΧΧΧ 


BEGIN READING THE DATA FROM A FILEDEF INTO ARRAYS. 


χχχχκχκχχχχκχκχκκκχχκχκχλκχκκχκχχκκαχκχχχχχκκχχχχχχκχχκχχκχκχχχχχχχκχκκχχκχχχκχχχ κκ 
ES 
IA 
READ(4,201) 55(Т), (RAT(I, μι ΡΕ ῬΡΕΝΕΆΝΤ(Σ,, 
AVE(I),TIME(I 
ООМО J=1, 10 


RATE(I, J)=RAT(I,J)-16 


100 CONTINUE 


S S T) E0.99) GO TO 50 
ου το 7 


50 CONTINUE 


大 大 大 大 大 类 大 类 大 大 大 大 类 大 炎炎 大 大 大 大 大 火 大 大 大 炎 大 大大 炎炎 大 大 大 大大 大 大 大 大 大 大 大大 大 类 大 大 大大 大大 大 大 炎 大 大火 大 大 大 大 大 大 大 


AT THIS POINT, THE DATA HAS ALL BEEN STORED IN THE ARRAYS SS(I), 
RATE(I,J),BASEL(I),CODE(I),ANT(I),RAVE(I), AND TIME(I). 


大 大 火 大 火炎 大 大 火 大 大 大 火 大 炎炎 大 大 大 大 大大 火炎 大 大 大 大 大 火 大 大大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 炎炎 大 大 大 大 大 火炎 大 大 大 大 


大 大 炎炎 类 大 大 大大 大 大 大 大 大 大 大 大 炎炎 大 类 大大 大火 大 大 炎炎 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 炎炎 大 大 炎炎 大大 炎 大 大 大 大 大 类 大 


DETERMINE LENGTH,BL, AND DATE AND BEGIN WRITING OUT THE DATA. 


ZA TT AT Te s 


DO 160 I=1,N 
ү с m. 
IF(BASEL 
IF(BASEL 
IF(BASEL 
IF(BASEL 
IF(BASEL 


.99) GO TO 500 
.EQ.1) GO TO 41 
"EOC2MCOTOSA2 
.Е0.3) GO TO 43 
.EQ.4) GO TO 44 


I4 I EA E EA CH 
[т] 
к 
ων 


.EQ.5) GO TO 45 
41 LENGTH(I)=1061.24 


BL=A 
DATE(I)=30 
“ο το 5 


42 LENGTH(I)=2417 .48 


BL=B 
DATE(I)=29 


GO TO 55 
43 LENGTH(I)-4083.10 


DES 
DATE(1)- 27 
GOTO 


44 TENen 7 = ОЛСО 


12245 


AAA 


OMMAaOa) 3M ОООО (О (О 


AAANANA A 


BL=D 
DATE(1)=28 
GO TO 55 

45 LENGTH(I)=9861.08 
BL=E 


HATE (1)=28 


ο το ο 
大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 六 大 类 大 大 本 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大 大 


BEGIN STATISTICAL CACCUCATTONGS: 
大 大 类 大 大 炎炎 大 大 大 大 大 类 大 类 大 类 大 类 大 大 火 大 炎炎 六 大 大 六 类 大 大 大 炎炎 大 大 大火 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大大 大 大 大 大 大 六 


55 RATSUM(I)=0.0 


50н(1)= 20.0 
DO 56 J=1,10 
а μμ т ост η) 
RERR(I,J)=RATE(I,J)-LENGTH(I) 
SUM(1)=SUM(I)+RERR(I,J) 
56 CONTINUE 
ee а 36 
ХВЕВВ(1)=5Ш0М(1)/10.0 
SUMM(I)=0.0 
De, SMO τα 
SUMM(I)=SUMM(1I)+((RATE(I,J)-XRATE(I))**2) 
57 CONTINU 


E 
VRATE (I= (SUMM(I))/9. 0 
SUMMM 
рО 58 151, 19 
SUMMM(I)SSUMMM(I)*((RERR(I,J)-XRERR(I))**2) 
58 CONTINUE 
VRERR(1)- (SUHMM(I))/9. 0 


59 ποτ) οὐ τουςς 
60 IF 5511. ПЕСО УО 
65 WRITE(7 202)BL, LENGTH(I), TDATE(I), TIMBCL) ,CODE(1) ,RAWE(T) 


大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 


TEST FOR THE ANTENNASTYPE TO BE USED IN THENBHBODER: 
A efc AR RR RR eh fh ho bj hh e he hehe 


IF 
IF(ANT(I).EQ.2) WRITE 
IF(ANT(I).EQ.3) WRITE 


WRITE (7 206) 
大 大 类 大大 类 关 大 大 类 大 大 大 大大 大 大 类 大业 大 大 大 大 大 类 大 类 大大 大 大 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


7,204 
205 





ANT 1) 59.2) WRITE 
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BEGIN ΜΕΙΤΙΝΘΕΙ ΤΗΕ ΤΑΝΕ 
大 大 大 大 类 类 大大 大 类 大 大 大 类 大大 大 类 六 大 大 大 六 大大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大 大 大 大 大 大 


80 CONTINUE 
| SS(I),XRATE(I),VRATE(I),XRERR(I) 
WRITE(7,209) RATSUM(I),sUM(I),sumM(ÍS,suMmMM(Í) 


160 CONTINUE 
大 类 大 大 大 大 六 大 大 大 类 类 大 大 大 大 六 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 大 大 


VRERR(I) 


FORMAT STATEHENTS: 
大 大 大 大 大 大 六 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


201 ER ДЕ 22 10/14 τσ, 11, 11 11 ἡ, 2, 

202 FORMAT(//2X, ! BASE LINÉ AL 2X, ‘LENGTH = | BIS METERS TE; 
1 “JANUARY АТ 1984',2Х,14 ÍX, HOURS' ,/, 18X, ' CODE: ТЕСЕУ 
1 'RANGE AVERAGE SETTING: μια 

203 FORMAT(22X,'ANTENNA TYPE: MEDIUM GAIN!) 

204 FORMAT(22X,'ANTENNA TYPE: OMNIDIRECTIONAL' ) 

205 FORMAT(22X,'ANTENNA TYPE: HIGH GAIN' ) 

206 FORMAT(/, 4X, 'SIGNAL' , 5X, |AVERAGH 6: eve RAGE SA; 
1 'RANGE ERROR',/,3X,'STRENGTH' ,6X,'RATE! ,6X, 'VARIANCE' , 5X, 
1 'RANGE ERROR! ,3X,'VARIANCE' ) 


126 


208 pe DA rs L AX Ph 1 SR F6 1) 
C 209 FORMAT(4X 'RATSUM= ',.F11.1,4X,'SUM- ',F6.1,4X,'SUMM- ',F15.1,4X, 
c 1 'SUHMM- ',F6.1) . 
500 CONTINUE 
STOP 


END 
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OOO OOO OOO OO OG а 


AAAN 


7 


50 


大 大 大 


大 大 大 


ΩΩΏΩΩΩΩ 


53 
60 
70 
80 


160 


大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


PROGRAM: SSPLOT 
AUTHOR: BRUCHESBCUHITIDBARD 
DATE: JUNE 6, 1984 
PROGRAM DESCRIPTION: C THIS FORDBANOCEBEERADSUSESCDESOSPERA 
SUBROUTINES TIO T CENERALE PLOTS OR 
MINI-RANGER SIGNAL STRENGTH VERSUS 
MEAN RANGE ERROR WITH STANDARD 
DEVIATIONS. 
INEUT: Ырыксыз FROM À DATA FILE REFERENCED WITH À 
EIE Be 
OUTPUT: THE RESULTANT PLOT CAN BE DIRECTED TO ONE OF THREE 
PLOTTERS: VERSATEC 8222A, VERSATEC O900A, OR A 
TEKTRONIX 4631 HARDWARE DEVICE. 


大 类 大 类 大 火 大 大 大 大 类 火炎 大 大 大 大 大 大 大 大大 类 大 大 火炎 大 大 大 大大 大 大 大 大大 大 大 大 大 大 大 大 大 大大 大 大 大 大 大 大 大 大 大 大 大 大大 大 大 大 大 大 大 大 大 大 


INTEGER I,J,K,L,M,SS,TIME,RAVE,CODE,ANT, BASEL, 

1 SSA,SS1,SS2,5S3,P,0,R,S,ANTA, IPAK 

DIMENSION 5 ANT 190) nn 500), 

1 SS(1500) , KRERR(1500), XRATE(1500) , VRERR(1500) , VRATE(1500), 

1 XRERR1(1500),XRERR2(1500),XRERR3(1500) ,XRERRÁ(1500), 

1 SSA(1500) ,S$1(1500) ,$S2(1500) ,SS3(1500) ,RAVE(1500),STD1A(1500), 
1 STD1B(1500) ,STD2A(1500) ,STD2B(1500) ,STD3A(1500) ,STD3B(1500), 

1 VRERRA(1500) , VRERR1 (1500) , VRERR2(1500) , VRERR3(1500) ,ANTA(1500), 
1 TEMP1(1500),TEMP2(1500), TEMP3(1500) , IPAK(100 

REAL XRERR,XRERRA,XRERRl,XRERR2,XRERR3,XRATE,VRATE,VRERR, 

1 STD1A,STDiB,STD2À,STD2B,STD3A,STD3B,VRERRA,VRERR1,VRERR2,VRERR3, 
1 TEMP1 ,TEMP2 ,TEIP3 


大 大 大 大 大 大 大 大 大 大 RERD IN THE STAT DATA FROM À FILEDEF. 


κ 
Bi 208) SS(I),XRATE(I),VRATE(I),XRERR(I),VRERR(I),BASEL(I), 
1CODE(I),ANT(I),RÀAVE(I),TIME(I) 

IF(SS(I1).EQ.99) GO TO 50 

GO TO 7 


CONTINUE 
ΝΞΙ 


大 大 大 大 大 大 大 AT THIS POINT, ALL THE DATA HAS BEEN STORED IN ARRAYS. 
AAKKKKK DETERMINE, BY COMPARISONS, THE DATA TO BE PLOTTED. 


DO 160 I=1,N 
IF(BASEL(I).EQ.5) GO TO 60 
GO TO 160 
IF(CODE(I).EQ.O) GO TO 70 
GO TO 160 
IF(RAVE(I).EQ.20)GO TO 80 
СО ТО 160 

SSA(P)=SS(I) 

ara AER (т) 

VRERRA(P)SVRERR(I 

ANTA(P)-ANT(I) 

P=P+1 


η πω H £ A > C. 
у» О О О ОО 


[j H H H H H H H H H H 


CONTINUE 
IG 


J=P-1 
ο ΑΝ. AT WWILSOPOINT ONLY IHE "BASE LINE A/CODE 0/RAVE=20"' 


GO GOO OOO 


ОООО ОО (О 


айда AAA AAA AAA 


OQ OO 


VALUES HAVE BEEN ACCEPTED INTO NEW ARRAYS 
πε ΠΙΤΕΡ BY ΤΗΕ INTEGER P. 


AX**AX***A4*** NOW LOAD THESE VALUES BY ANTENNA-TYPE INTO THREE 
SEPARATE ARRAYS AND CALCULATE THEIR STANDARD 
DEN MOTIONS ERON TEEIR MEANS. 


DO 180 P=1,J 








IF(ANTA(P).EQ.1) GO TO 191 
IF(ANTA(P).EQ.2) GO TO 192 
IF(ANTA(P).EQ.3) GO TO 193 
191 551(0)-55А(Р) 
Е 
VRERR1(0)=VRERRA(P 
STDIA()“HRERRI (0)+s | 
Т О)=ХВЕВВ1(О)-50ЕТ(УВЕВВ1 (0 
=0+ 
GO IO 180 
192 SS2(R)=SSA(P) 
s. ) “sine, 
VRERR2(R)=VRERRA(P 
r М МЕЗ ss ОВ 人 
STD2B(R)=XRERR2(R)-SORT(VRERR2(R 
R=R+1 
GO TO 180 
193 SS3(S)=SSA(P) 
ο ο πε 
VRERR3(S)=VRERRA(P 
2009262 S м цы БИ 
STD3B(S)=XRERR3(S)-SQRT(VRERR3(S 
S=S+1 
GO TO 180 
180 CONTINUE 
К=0-1 
L-R-1 
M=S-1 


大 大 六 大 大 大 六 大 大 大 大 大 大 PERFORM THE NECESSARY PLOT/SETUP ROUTINES. 
六 大 大 大 大 大 大 大 大 大 TINITIALIZE PLOTTER AND SET PAGE SIZE. 
CALL PRTPLT(72,6) 
CALL COMPRS 
CALL VRSTEC (0,0,0) 
CALL TEK618 
CALL PAGE (11.,8.5) 
CALL BLOWUP(1.25) 
大 大 大 大 大 大 大 大 大 大 MAKE SS AXIS IN INTEGERS 
CALL XINTAX 
大 大 大 大 大 大 大 大 大 大 ORIENT Y AXIS NUMBERS HORIZONTALLY. 
CALL YAXANG (0.0) 
AX******** SELECT TOTAL PLOT AREA. 
CALL ARÉA2D(8.,6.) 
大 大 大 大 六 太太 大 大 大 SCALE AND LABEL AXES WITH SELF-COUNTING OPTION. 
CALL YNAME('RANGE ERROR (M) P0 
CALL XNAME се STRENGTHS' ,100 
CALL BASALF('STAND') 


大 大 大 大 大 大 大 大大 大 CREATE THE HEADING AND SET THE PRINT TYPE. 


129 


ООО О OQ O O ООО 


ОО О О О О О О О О О О О 


QQ 


СУСУСУ 


CALL SWISSM 
CALL HEADIN ('MINI-RANGER III BASELINE CALIBRATIONS',100,1.5,1) 


大 大 大 大 大 大 大 大 大 大 DEFINE BOTH AXES AND FRAME THE SUBPLOT AREA. 


CALL DUPLX 

CALL GRAF(18.0 ,1-0,3.0,-10.0,2.0,10.0) 
CALL THKFRM(.02) 

CALL FRAME 


大 大 大 大 大 大 大 大 大 大 SET UP LEGEND DATA AND TEXT FOR LEGEND NAME. 


CALL LINESP (2.0) 

CALL LINES ('MEDIUMS' , IPAK,1) 
CALL LINES OUNTO TEAK, 2) 
CALL LINES ('HIGHS''IPAK,3 


大 大 大 大 大 大 大 大 大 大 FIND SIZE FOR LEGEND. 


ҮН=ҮГЕӨМО (ТРАК, 3 
YW=YLEGND ( IPAK,3 
CALL MYLEGN(' ANTENNA TYPE:$' PLC) 
CALL LEGLIN 
AKKKKKKKKK DEFINE CURVE TYPES AND EXECUTE THE CURVES. 


kkkkkkkkkk PLOT THE MEDIUM GAIN ANTENNAS TURSE 


CALL MARKER(15) 
CALL DOR 
CALL THKCRV(. 025) 
DO 195 Dni. 
ТЁИР1 (о)= FLOAT(SS1(Q)) 
195 CONTINUE 
CALL CURVE(TEMP1,XRERR1,K,1) 


大 大 大 大 大 大 大 大 大 大 大 PLOT THE OMNIDIRECTIONAL ANTENNA VALUES. 


CALL MARKER(16) 
CALL DASH 
CALL THKCRV(.025) 
DO 196 R=1,L 
TEMP 2(R)=FLOAT(SS2(R)) 
196 CONTINUE 
CALL CURVE(TEMP2,XRERR2,L,1) 


大 大 大 大 大 大 大 大 大 大 大 大 PLOT THE HIGH GAIN ANTENNA VALUES. 


CALL MARKER(17) 
CALL RESET('DASH' ) 
CALL THKCRV(.025) 
DO 197 S=1,M 
TEMP3(S)=FLOAT(SS3(S)) 
197 CONTINUE 
CALL CURVE(TEMP3 , XRERR3 ,M,1) 


大 大 大 大 大 大 大 大 大 大 OTHER PLOT/LEGEND INFORMATION: 


CALL MESSAG( ‘BASELINE BS! .100 3.506055 -- 

CALL MESSAG('RANGE AVERAGE: 20$!',100,3.50,5.25 

CALL MESSAG('RANGE AVERAGE: "4095. 1007 998805925 

CALL MESSAG('RANGE AVERAGE: 75$',100,3.50,5.25 

CALL MESSAGC CODESNOS ον Ἱοο o n 

CALL MESSAG ТЕБЕН ΝΟ, 15 007575 589 

CALL MESSAG('NOTE: VERTICAL BÁRS ÁRE STANDARD DEVIATIONS. SE 
l ТОО 5 50) 


大 类 大 大 大 大大 大 类 大 DRAW THE LEGEND. 
CALL LEGEND(IPAK,3,1.0,3.75) 
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ОООО 


AQAA ANQ 


NQ 


эсен PLOT VERTICAL BARS (STANDARD DEVIATIONS). 
CALL A EM 
CALL BARPAT(16) 
CABESVBARSCTEMPI,SIDIB,STDIA,K 
CABEB-VBARS(TEHMP2 STDZBPSIDZA.L 
CALE VBARSUTENUPSTSTDSB,SIDSATM 
大 大 大 本 大 大 大 大 x READ FORMAT STATEMENT. 
208 ОЛ 12/05 ЕБІ το ος δι 1LLXRRS .1 ,3X,IL,I1,I1,I2,I4) 
W********* CLOSE OUT METAFILE AND RETURN TO LEVEL ZERO. 


CALL ENDPL(0) 
CALL DONEPL 





SLOR 
END 
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大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


PROGRAM: WALK 
AUTHOR: BRUCE F. HILLARD 
DATE: AUGUST 24, 1984 
PROGRAM DESCRIPTION:  IHISSEORIEADIMPROGBADNUSESSDISSPEA 
SUBROUTINES TO GENBEAIE PEOTS” OF 
BASE LINE DISTANCE VERSUS 
MEAN RANGE ERROR. SIGNAL STRENGTH 
VALUES ARE USED AS MARKERS. 
INPUT: UST CONES FROM A DATA FILE REFERENCEDIWITH A 
ΤΡ. 
ОСТРОЕ: THE RESULTANT PLOT CAN BE-DIRECTED TO ONE OF THREE 
PEOTTERS: VERSATEC 2 ορ ο 088, ο À 
TEKTRONIX 4631 HARDWARE DEVICE. 


大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


AAAAAAAAAAAAAAAAAAANA 


INTEGER I,J,K,TIME,RAVE,CODE,ANT,BASEL,SSA,SS,BASELA,B 

1N,NN,NNN, NNNN,M,MM, MMM,MMMM, P, PP,PPP,PPPP,R,RR,RRR 

DIMENSION το. ο, 

1 SS(1500), XRERR(1500),XRATE(1500) , VRERR(1500) , VRATE(1500), 

1 XRERRA(1500),SSA(1500) , RAVE(1500) , BASELA(1500) 

REAL XRERR, XRERRA, XRATE, VRATE, VRERR, L 

1 Х4,ХХ4,Х5,ХХ5,Х6,ХХ6,Х7,ХХ7,Х8,ХХ8,Х9,ХХ9,Х10,ХХ10, 

! X11, XXÍ1,X12,XX12,X13,XX13,X14,XK14,XÍ5,XX15,%16,X%16,X17,XX17, 
A 


C 
C ***X******* INITIALIZE PLOTTER AND SET PAGE SIZE. 
C 
5 CALL PRTPLT(72,6) 
CALL COMPRS 
е CALL VRSTEC (0,0,0) 
C CALL TEK618 
CALL PAGE (11.,8.5) 
© 
С 
C 大 大 大 大 大 大 大 大 大 大 ORIENT Y AXIS NUMBERS HORIZONTALLY. 
С 
CALL YAXANG (0.0) 
Ç 
C ********** SELECT TOTAL PLOT AREA. 
C 
CALL AREA2D(8.,6.) 
С 
C 大 大 大 大大 大 大 类 大 大 SCALE AND LABEL AXES WITH SELF-COUNTING OPTION. 
C 
CALL YNAME('RANGE ERROR (M) ion 
CALL XNAME Е LINE DISTANCE (M)$',100) 
CALL BASALF('STAND' ) 
C 
C W********* CREATE THE HEADING AND SET THE PRINT TYPE. 
С 
CALL HEADIN ('BASE LINE DISTANCE VS. RANGE ERRORS',100,1.5,1) 
C 
C ********** DEFINE BOTH AXES AND FRAME THE SUBPLOT AREA. 
C 
CALL DUPLX 
CALL СБАЕ(0,0 10900,0 102300900 
CALL THKFRM(.02) 
CALL FRAME 
С 
C 
C 大 大 大 大 大 大 大 大 大 大 OTHER PLOT/LEGEND INFORMATION: 
С 


CALL ἂν ο ТОНЕ πο, ОП ЕШ 850, 252225) 
CALL MESSAG('NUMBERS ARE SIGNAL STRENGTHS $', 
1 100 2250 5550) 


152 


С 
e 
C A**WA^ANAA* READ IN THE STAT DATA FROM A FILEDEF. 


C 
J=0 
1-0 
7 I=I+1 
p s SS(I),XRATE(I),VRATE(I),XRERR(I),VRERR(I),BASEL(I), 
1CODE(I) ANT(I) ,RAVE(I), TIME(I) 
IF (SS (1). ΕΟ.994 σο ΤΟ 50 
GO TO 
50 CONTINUE 
N-I 
c 
ς 大 大 大 大 大 大 大 大 大 大 RERAD FORMAT STATEMENT. 
S00) Posten ος 12 ?Z Б6.1,8Х,Е3.1,10Х,Ғ4.1,11Х,Е3.1,3Х,11,11,11,12,14) 
ς 大 大 大 大 大 大 大 大 大 大 AT THIS POINT, ALL THE DATA HAS BEEN STORED IN ARRAYS. 
С KXkXkkX***** SORT OUT THE DESIRED CODE. 
C 
ро 400 І=1, 
tÉ (CODE (1). .NE.0) GO TO 400 


SSA(J)=SS(I) 
ОЕ 


HH 


BASELA (J )=BASEL ) 
400 CONTINUE 


C 
ο KK IS THE TOTAL NUMBER OF POINTS TO BE PLOTTED. 


C 
K=J 
B=0 


5 
C 大 大 大 大 大 大 大 大 大 大 大 大 大 大 DO LOOP THROUGH THE DATA ONCE PER BASE LINE. 
500 CONTINUE 
B=B+1 
Е IF(B.EQ.6) GO TO 2000 
C 大 大 大 大 大 大 大 大 大 大 大 大 大 大 SET ALL COUNTERS TO ZERO. 


C 
N=0 
A4-0.0 
XX4=0 .0 
NN=0 
Х5-0.0 
XX5=0.0 
NNN=0 
X6=0.0 
XX6=0 .0 
NNNN=0 
X7=0.0 
XX7=0.0 
Iz 
X8-0.0 
XX8=0.0 
MM=0 
Х9-0.0 
XX9=0 .0 
MMM=0 
Х10-0.0 
XX10=0 .0 
MMMM=0 
X11=0.0 
ХХ11=0.0 
Р=0 
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Х1250.0 
XX12=0.0 
PP=0 
X13=0.0 
XX13=0.0 
PPP=0 
Χ14Ξ0.0 
XZ14=0.0 
PPPP=0 
X15=0.0 
АХ15=0.0 
R=0 
X16=0.0 
XX16=0.0 
RR=0 
X17=0.0 
XX17=0.0 
RRR=0 
X18=0 .0 
ХХ18=0.0 


C 
C 关 大 大 大 大 大 大 大 大 大 大 大 大 大 SORT THE DATA BY BASE LINES AND SIGNAL STRENGTHS. 


e 


DO 1000 J-1,K 


С 
C 大 大 大 大 大 大 大 大 大 大 大 大 大 大 SORT FOR ONE BASE LINE AT A TIME. 


IF(BASELA(J).NE.B) GO TO 1000 


С 
CAARAAKAARAKAAK SORT THE SIGNAL STRENGTHS. 


C 
100 
TG 
102 
105 
104 
105 
106 


107 


CONTINUE 

HE απ IOl 
X4=X4+XRERRA (J 

N=N+1 

GO TO 1000 

CONTINUE 

Ἢ, Ἡ СО ТО 102 
X5=X5+XRERRA(J 

NN=NN+1 

GO TO 1000 

CONTINUE 

| πο πο 5 
X6=X6+XRERRA(J 

NNN=NNN+1 

GO TO 1000 

CONTINUE 

E GO TO 104 
X72X7*XRERRA(J 
NNNN=NNNN+1 

GO TO 1000 

CONTINUE 

rU CO TO 105 
X8=X8+XRÉRRA(J 

M=M+1 

GO TO 1000 

CONTINUE 

PNE GO TO 106 
X9=X9+KRERRA(J 

MM=MM+1 

GO TO 000 

CONTINUE 

DE (SSH ΙΕ 10) GO TO 107 
X10=X10+ZRERRA(J) 
MMM=MMM+1 

GO TO 1000 

CONTINUE 

IF(SSA(J).NE.11) GO TO 108 
X11=X11+XRERRA (J) 
MMMM=MMMM+1 
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GO TO 1000 
108 CONTINUE 
IF(SSA(J).NE.12) GO TO 109 
X12=X12+XRERRA(J) 
P=P+1 
GO TO 1000 
109 CONTINUE 
IF(SSA(J).NE.13) GO TO 110 
X13=X13+XRERRA (J) 
PP=PP+1 
GO TO 1000 
110 CONTINUE 
IF(SSA(J).NE.14) GO TO 111 
X14=X14+XRERRA (J) 
PPP=PPP+1 
GO TO 1000 
111 CONTINUE 
SSA(J).NE.15) GO TO 112 
-X15*XRERRA(J) 
Pppp PPPP+1 
GO TO 1000 
112 CONTINUE 
IF(SSA(J).NE.16) GO TO 113 
XX16=X16+XRERRA(J) 
R=R+1 
GO TO 1000 
113 CONTINUE 
IF(SSA(J).NE.17) GO TO 114 
X172X174XRERRA(J) 
RR=RR+1 
GO TO 1000 
114 CONTINUE 
IF(SSA(J).NE.18) GO TO 1000 
к= =K18+XRERRA (J) 
RRR=RRR+1 
1000 CONTINUE 


C 
C*A*X*X*X**AX«4aX COMPUTE THE. MEAN. RANGE ERROR FOR EACH SIGNAL STRENGTH 
C PUD ASSIGN) THE BASE LINE LENGTH PO ΡΕ PLOTTED: 


C 
IF(B.NE.1) GO TO 90 
L=1061.24 
GO TO 600 
90 IF(B.NE.2) GO TO 91 
L=2417.48 
GO TO 600 
91 IF(B.NE.3) GO TO 92 
L=4083.10 
GO TO 600 
92 IF(B.NE.4) GO TO 93 
EFIE? 
GO TO 600 
93 CONTINUE 
L-9861.08 
600 CONTINUE 
Ckkkkkkkkkkkkkk COMPUTE THE MEAN RANGE ERRORS. 


С 
XX4-X4/FLOAT (N) 
XX52X5/FLOAT(NN) 
kX6=X6 / FLOAT (NNN) 
XX7=XK7 / FLOAT (NNNN) 
XX8=K8 /FLOAT(M) 

XX9=X9 / FLOAT (MM) 
XX10=X10/FLOAT (HMM) 
XX11=X11/FLOAT (MMMM) 

ХХ12-Х12/ЕІОАТ(Р) 
XX13=X13/FLOAT(PP) 
XX14=K14/FLOAT(PPP) 
XX152X15/FLOAT(PPPP) 
XX16=X16/FLOAT(R) 


B35 


UN 
XX18-2X18/FLOAT(RRR) 


C 
CRARKKAKKARAKAK SET THE CHARACTER HEIGHT. 
CALL HEIGHPO O ПС), 
C 
C 大 大 大 大 大 大 大 大 大 大 大 大 大 大 PLOT THE VALUES FOR ONE BASE LINE. 


Ç 

ТЕ(ХХ4.ЕОШО 0 БСО πο 201 
CALL RLINT(4,L,XX4) 

201 CONTINUE 
IF(XX5.EQ.0.0) GO TO 202 
CALL RLINT(5,L,XX5) 

202 CONTINUE 
IF(XX6.EQ.0.0) GO TO 203 
CALL RLINT(6,L,XX6) 

203 CONTINUE 
IF(XX7.EQ.0.0) GO TO 204 
CALL Επ Л К 5 

204 CONTINUE 
IF(XX8.EQ.0.0) GO TO 205 
CALL REINT(B, L; XX8) 

205 CONTINUE 
IF(XX9.EQ.0.0) GO TO 206 
CALL RULIN ООО -Ὁ 

206 CONTINUE 
IF(XX10.EQ.0.0) GO TO 207 
CALL RLINT(10,L,XX10) 

207 CONTINUE 
ИБО ὁ ο). σϱ.το. 205 
CALL REINT(11,L, #11) 

208 CONTINUE 
ТЕТІКТІ) CO TO 209 
CALL RLINT(12,L,XX12) 

209 CONTINUE 
TF (30013 E90, 0) GO ТӘ 210 
CALL RLINT(13,L,XX13) 

210 CONTINUE 
J (O 12069.0720) GO TO 241 
C PP BETNT( IAE XI) 

211 CONTINUE 
IF(XX15.EQ.0.0) GO TO 212 
CALL RLINT(15,L,XX15) 

212 CONTINUE 
IF(XX16.EQ.0.0) GO TO 213 
CALL RLINT(16,L,XX16) 

213 CONTINUE 
Pex 7eneno.0) Со TO 214 
Cink LT 1? 1. ХХ1Т) 

214 CONTINUE 
IF(XX18.EQ.0.0) GO TO 215 
CALL RLINT(18,L,XX18) 

215 CONTINUE 

C 


GO TO 500 
2000 CONTINUE 


C 
CAAAAAAKAAAAKA CLOSE OUT METAFILE AND RETURN TO LEVEL ZERO. 
C 

CALL ENDPL(0) 

CALL DONEPL 


STOP 
END 
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ЖАХААХАХХАААААХАААХАХАКААААХАХХАХАХААХАХАХААААХАХАААХААХАХАХАХАХАХАХАААААХАХАЖ 


PROGR: POEYEIT 

AUTHOR: BRUCE F. HILLARD 

DATE: SEPTEMBER 1984 

PROGRAM DESCRIPTION: THIS FORTRAN PROGRAM, TAKEN FROM 
GERALD'S NUMERICAL ANALYSIS BOOK, 
Εις ΕΕ ΝΟΚ ουν UP TO DEGREE 
7 TO THE MEAN RANGE ERROR VALUES 
FORIMINI -RANGER CODES. 


INPUT: MEAN RANGE ERRORS ARE INCLUDED IN A DATA STATEMENT. 
OUTPUT: THE RESULTS ARE WRITIEN TO A PRINTER OR DISK FILE. 


大 大 六 火炎 大 大 大 大 大 大 类 大 大 大 大 大 类 大 大 大 大 大 火 大 大 火炎 火炎 交大 大 大 大 大 大 大 大 大 大 火 火 火 火 大 六 六 大 大 火 火 大 大 大 大 大 火 太 大 火 六 大 大 大 大 大 火 大 


AAAAAAAAAAAAAAAANANAA 


IMPLICIT REAL*4 (A-H,0-Z) 
DIMENSION X(100),Y(1 60). C(100), (10.1 


nud 
1:06) 20009 mS CO) T RQ9)TX 00) , X (119, 


DATA XO), 0297 x(3), X(4), 2555) 
1 Тм 
IE) О ο τς ο το 13.0 12 0 1150 10.0/920,8.0,7.0, 
] 6.0,5.0,4.0/ 
р DATA 119,70); ООЛ πα το ο ποτ ο στ) πο) (т, 
ү ,Y | 
1 n 4,1.4,0.8,-0.3,0.5,0.6,0.8,0.5,0.4,0.0,-0.5,-0.9,0.5, 
575! 
c 
C THIS PROGRAM IS USED TO IN FITTING A POLYNOMIAL TO A SET OF DATA. 
C THE PROGRAM READS IN N PAIRS OF X AND Y VALUES AND COMPUTES 
C THE COEFFICIENTS OF THE NORMAL EQUATIONS FOR THE LEAST SQUARES 
C METHOD. 
C PARAMETERS ARE: 
5 ΤΥ ARRAY ΟΕ X AND Y VALUES. 
C N NUMBER OF DATA PAIRS 
C MS ,MF THE RANGE OF DEGREE OF POLYNOMIALS TO BE COMPUTED. 
С THE HAXIMUH DEGREE IS 9. 
С А AUGMENTED ARRAY OF THE COEFFICIENTS OF THE NORMAL 
С EQUATIONS. 
c Ç ARRAY OF COEFFICIENTS OF THE LEAST SQUARES POLYNOHIALS. 
С 
N=14 
MS=1 
HF=7 
C COMPUTE MATRIX OF COEFFICIENTS AND R.H.S. FOR HF DEGREE. 
C HOWEVER, FIRST CHECK TO SEE IF MAX DEGREE REQUESTED IS TOO LARGE. 
C IT CANNOT EXCEED N-1. IF IT DOES, REDUCE TO EQUAL 
C N-1 AND PRINT MESSAGE. | 


BECOME.LE.(N-1)) GO TO 5 
MF=N-1 
WRITE(6,200) MF 
200 FORMAT (/, 'DEGREE OF POLYNOMIAL&CANNOT EXCEED N-1',//, 
| x UESTED MAXIMUM DEGREE OF POLYNOMIAL TOO LARGE-REDUCED πο. 
5 ΗΕΕΙΞΗΕΓΙ 
МЕР2=МЕ+2 
C PUT ONES INTO A NEW ARRAY. THIS WILL HOLD POWERS OF THE X VALUES 
C AS WETEROCEED. 
DO 10 I=1,N 
10 AN(I)=1 
C COMPUTE FIRST COLUMN AND N*1 ST COLUMN OF A. I MOVES DOWN THE 
Ου US OVER THE N VALUES. 
DOSSOSI-I/HMEFPI 
ҮЛГІ. 
А(І,МЕР2)-0.0 
DO 20 J=1,N 


1” 7 


А ELE а(1, 1) + XN(J) 

А(І,МЕР2)- А(І, МЕР2) Fo 

20 XN(J)= XN( J) (9) 

30 CONTINUE 

C COMPUTE THE LAST ROW OF A. I MOVES ACROSS THE COLUMNS, 
C J SUMS OVER N VALUES. 


(в 
DO 50 I=2,MFP1 
A(MFP1,I) = 0. 0 
DO 40 Ј = 1. 
40 XN(J) = RN(S) * X(J) 
50 CONTINUE 
С МОЙ ЕТЕР ІМ ЕНЕ БЕЗБЕН ТИЕ ГА МАТЕ: I MOVES DOWN ЕНЕ 
C ROWS, J MOVES ACROSS THE COLUMNS. 
DO 70 ξ 2,MFFP1 
DO 60 І = 1 ΜΕ 
60 A(I,J) = A(I*1,J-1) 
70 CONTINUE 
C WRITE OUT THE MATRIX OF NORMAL EQUATIONS. 
WRITE(6,201) (CA (1,J), J = 1,MEP2), I - І ШЕБІ 
201 FORMAT(/, κι» Е 
C NOW CALL A SUBROUTINE TO SOLVE THE SYSTEM] =0C THIS ЕОҚ ЕДӨРИЗЕСЕНРЕ 
C FROMSMS ο n 
С СЕТ THE ο Ομ ποστ ОМ OF A 
CALL LUDCMQ (A, ΜΕΡΙ, 10) 
C RESET THE R.H.5. INTO C. WE NEED IO DOS PHPESSPORVESGH DECREE. 
HSPI1 =S + 1 
DO 95 I = MSP1,MFP1 
B18) S10) eS ily r 
90 C(J) = A(J,MFP2) 
CALL SOLNQ (A, C, I, 10) 
THI = T 1 
C NON WRIIE OUT THE μος, 955 D. LEAST SQ POLYNOMIAL. 
WRITE(6,202) IM1, (C(J), uj 
202  FORMAT(//, "κα, 'FOR DEGREE oF М) COEFFICIENTS ARE: ',/,10X, 
C COMPUTE AND КЕ THE VALUE OF BETA = SUM OF DEVIATIONS SQUARED, 
C WHICH IS = (N - 1 - M). 
BETA=0.0 
DO 94 IPT = 1,N 
SUM=0 .0 
DORIS ICOBE = >P I 
JCOEF- I-ICOEF + 2 
SUM=(SUM + C(JCOEF))*X(IPT) 
93 CONTINUE 
SUM=SUM+C(1) 
BETA=BETA + (Y(IPT) - SUM)**2 
94 CONTINUE 
BETA=BETA/(N-1) 
WRITE(6,203) BETA 
203 FORMAT(1H1,10X, ' BETA IS ',F10.5) 
95 CONTINUE š 
STOP 
END 


SUBROUTINE LUDCMO (A, N, NDIM) 
DIMENSION A(NDIM,NDIM) 
C THIS SUBROUTINE FORMS THE LU EQUIVALENT OF THE SQUARE 
C COEFFICIENT MATRIX A. THE LU, IN COMPACT FORM, IS RETURNED 
C IN THE A MATRIX SPACE. THE UPPER TRIANGULAR MATRIX U HAS ONES 
C ON ITS DIAGONAL - THESE VALUES ARE NOT INCLUDED IN THE RESULT. 
DO 30 ТТ 
DO 30 J =2,N 
SUM=0. 
IF(J.GT.1). GO TO 15 
ΠΠ = 
рО 10 К=1, T 
10 SUM = SUM + A(T, K)AA(K, Т 
A(1. J) = А(Т ee od 


w38 


GO πο 30 
15 IM] = I 
IF(IM1. EQ. ὃ) ιο 29 
DO 20 K 5 1,IMl 
20 SUM = SUM 'A(I,K)%*A(K,J) 
C TEST FOR SBALL VALUE ON DIAGONAL. 
25 ТЕ(АВ5(А(Т, 1)); LT.0.0000001) GO TO 99 
А(І,Ј) = (А(І,Ј) - SUM)/A(I,I) 
30 CONTINUE 
RETURN 
99 WRITE(6,100) I 
100 FORMAT(/,' REDUCTION NOT COMPLETED BECAUSE 
1 SMALL VALUE FOUND FOR DIVISOR IN ROW ',I3) 
RETURN 
N 


SUBROUTINE SOLNQ (A,B,N,NDIM) 
DIMENSION A(NDIM, NDIM), B(NDIM) 


C THIS SUBROUTINE FINDS THE SOLUTION TO A SET OF N LINEAR 
C EQUATIONS THAT CORRESPONDS TO THE RIGHT HAND SIDE VECTOR B. 
C THE A MATRIX IS THE LU DECOMPOSITION EQUIVALENT TO 
C THE COEFFICIENT MATRIX OF THE ORIGINAL EQUATIONS, AS 
C PRODUCED BY LUDCMP. THE SOLUTION VECTOR IS 
C RETURNED IN THE B VECTOR. 
C DO THE REDUCTION STEP 
B) "BOE AGI) 
DO 20 I 2 2.N 
ΜΙ =I- 1 
SUM = 0.0 
DO 10 K = 1,IM1 
10 SUM = SUM * A(I,K)^B(K) 
20 B(I) = (B(I) - SUM)/A(I,I) 
C NOW WE ARE READY FOR BACK SUBSTITUTION. REMEBER THAT THE 
C ELEMENTS OF U ON THE DIAGONAL ARE ALL ONES. 
DO 40 J = 2,N 
NHJP2 = N - J +2 
NHJP12N-J-*1 
SUM = 0. 
DO 30 K = NMJP2,N 
30 SUM = SUM + А(ШЫР1, K)*B(K) 
40 B(NHJP1) = B(NMJP1) - SUM 
RETURN 


END 


5» 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkĘk 


PROGRAII: ROLY 
AUTHOR: BRUCE F. HILLARD 
DATE: SEPTEMBER 1984 
PROGRAM DESCRIPTION: THIS EORTRANCEROGRAMQGUSESODISSPEA 
OSUBROUTINESSEO CENERATE PECTS OF 
MINI-RANGER RANGE ERRORS THAT ARE 
PREDICTED FROM POLYNOMEAL REGRESSION 
CURVES. ALSO РГОТТЕ Б НЕ THE 
MEAN RANGES ERRORS USED TO CREATE 
THE REGRESSION CURVES. THE RANGE 
ERRORS ARE ASSUMED TO BE A FONCTION 
ΟΡ ΤΟΝΟΙ ολο η, 
ΙΝΕ: PLOT VALUES ARE GENERATEDOWIIIHIINOLHE PROGRAMSBY 
IHE REGRESSION CURVES. ІНЕ ΠΕΑΝ RANGE ERRORS ARE 
INPUT THROUGH A DATA STATEMENT. 
OUTPUT: THE RESULTANT PLOT CAN BE DM IED TO ONE OF THREE 
PLOTIERS® VERSATEC S222 ες τις οσο ο 5 
TEKTRONIX 4631 HARDWARE DEVICE. 


大 大 大 大 大 大 大 大 大 大 大 大 类 大 大 大 类 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 大 


СУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУ 


IMPLICIT REAL*4 (А-Н,0-2) 
DIMENSION X(14),Y(14], Ү1(100) 72(100), ¥3(100) ,X1(100), 
00) ,X3(100) , IPAK(100 
),Х(4),Х(5),Х(6),Х(7),х(8),Х(9),Х(10),Х(11), 
,7.0,8,50,9202Л0011 000140000 
,Y(4),Y(5), (БУУЛ (2 τι ο 0 hh 1) 


;,71.4,-1.3,-22086,-1995090 0 m ο ο το... 


- ж 


大 大 大 大 大 大 大 大 大 大 大 大 大 PERFORM THE NECESSARY PLOT/SETUP ROUTINES. 
大 大 大 大 大 大 大 大 大 大 INITIALIZE PLOTTER AND SET PAGE SIZE. 
CALL PRTPLT(72,6 
CALL COMPRS 
CALL VRSTEC (0,0,0) 
CALL TEK618 
CALL PAGE (11.,8.5) 
CALL BLOWUP(1.25) 
大 大 大 大 大 大 大 大 大 大 MAKE SS AXIS IN INTEGERS 
CALL XINTAX 
大 大 大 大 大 大 大 大 大 大 ORIENT Y AXIS NUMBERS HORIZONTALLY. 
CALL YAXANG (0.0) 
大 大 大 大 大 大大 大 大大 SELECT TOTAL PLOT AREA. 
CALL AREA2D(8.,6.) 
Ak*k*k*k*k* SCALE AND LABEL AXES WITH SELF-COUNTING OPTION. 
CALL YNAME('RANGE ERROR (M) DO 


CALL XNAME E STRENGTH$' ,100 
CALL BASALF('STAND!) 


AQAA AAA AAA AAAA AA Οσο ο ο 


大 六 大 大 六 六 六 六 大 六 DERINE BOTH AXES AND FRAME THE ΡΜ. 
CALL DUPLEX 


СУСУСУСУ 
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(ОПП СЕДЕ( 18.0 анайы дс; 
САРЫ ІНЕБЕМ(С:02) 


CALL FRAME 
C Е 
в С CREATE THE LEGEND THAT IDENTIFIES ЕАСН PLOT. 
С 
С 
c kkkkkkkkkk OTHER PLOT/LEGEND INFORMATION: 
С 
С 
C**** COMPUTE THE FIRST DEGREE VALUES. 
C 
9 А1=-0.883 
B1=0.079 
ТЕМР1=4.0 
20 CONTINUE 
DO 200 I=1,65 
q IP1 
Yl(I Alt (BI*X1(1)) 
ТЕМР1=Х1(1)+0.2 
200 CONTINUE 
c 
C**** COMPUTE THE SECOND DEGREE VALUES. 
C 
A2=4.540 
B2=-1.133 
C2=0.058 
С 
ТЕМР2=4.0 
30 CONTINUE 
DO 300 J=1,65 
gu EM 
Y2 -h2 (82432 (3) )e (CZ^2(3) ^92) 
TEMP2=X2(J)+0 .2 
300 CONTINUE 
C 
C 
C**** COMPUTE THE THIRD DEGREE VALUES. 
C 
A3217.337 
B3=-5.661 
C3=0.530 
А D3=-0.015 
ТЕМР3=4.0 
40 CONTINUE 
DO 400 Kel 65 
KOS 
Y3 KSAS (oem 0 )«(c2*12 (6) 2) (n3^13 (0) 443) 
TEMP 3=X3(K)+0.20 
400 CONTINUE 
С 
С 大 大 大 大 大 大 大 大 大 大 大 SET UP LEGEND INFORMATION. 
С 
MAXLIN=LINEST(IPAK,100,25) 
1 CALL LINESP(2.0 
2 CALL LINES('COMPUTED MEAN VALUESS' , IPAK,1) 
3 CALL LINES('FIRST DEGREE VALUESS' , IPAK,2) 
4 CALL LINES('SECOND DEGREE VALUESS!,IPAK,3) 
5 CALL LINES('THIRD DEGREE VALUESS' , 1PAK,4) 
6 CALL MYLEGN('CURVE TYPES:',12) 
7 CALL LEGLIN 
C 
C 
C**** PLOT THE MEAN ERRORS AND THE COMPUTED VALUES. 
C 
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ΩΩ 


ΩΩωΩΩΩΩΟΩΩΩΩΩΩΩΩΩΩΩΩΩΩ ΩΩΩΩΩΩΩΩΩ а 


ΩΩ 


CALL THKCRV (029) 
CALL CURVE(X,Y,14,0) 


CALL THKERV( 0a 
CALL DOT 
CALE CURVE( A], Vie comer 


CALL THKCRV(.025) 
CALL DASH 
CALL CURVE(X2,Y2,65,0) 


CALL THKCRV(.025) 
CALL CHNDOT 
CALL CURVE(X3,Y3,65,0) 


CALL LEGEND(IPAK,4,1.0,4.0) 


RARAAKAAKK ALL CURVES HAVE BEEN PLOTTED. 
AX*X******* NOW COMPUTE THE M.A.S.S. AND PLOT IT. 


CLASTSX3 (2) 
DO 500 K=2,65 
CHECK=ABS (X3(K)-CLAST) 
IF(CHECK.GE.2.0) GO TO 600 
CLAST=X3 (K) 
500 CONTINUE 


хх ек Μο αι ΠΕ Ηλ ο 


600 CONTINUE 
D 
ХЕ-Х3(1.) 
ΧΤΞΧΕ 
E rc 
YT-Y3(L)40.1 


CALL VECTOR(XF,YF,XT,YT,1303) 
ХХ-ХЕ-0.3 
YY=YF40 Z 
САСО ΚΕΠ 
CALL MESES 
CALL MESSAG 


ο ος 1ου, τ, 
(CODE NO. το осе 
"CODE NO. 19',T007 50405590 





kkkkkkkkkk CLOSE OUT METAFILE AND RETURN TO LEVEL ZERO. 
CALL ENDPL(0) 
CALL DONEPL 
STOP 
DEBUG TRACE 
ΠΠ 


TRACE ON 
END 
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大 类 类 大 大 大 大 火 火 类 大 大 大 火 火 类 大 大 火 火 火 火 火炎 火炎 火 火 火炎 大 大火 大 类 火 大 大 大 类 类 火炎 大 六 炎炎 大 大 六 大火 大 火炎 类 火炎 六 六 类 火炎 类 六 大 火 六 大 


PROGRAM: BASELINE 
AUTHOR: BRUCE PF. HELEARD 
DATE : AUGUST 1985 
PROGRAM DESCRIPTION: THIS FORTRAN PROGRAM USES DISSPLA 
SUBROUTINES TO GENERATIE PLOTS OF 
MINI-RANGER RANGE ERROR VERSUS 
SIGNAL STRENGTH BY BASE LINE. 
ΠΟΙ: COMES FROM À DATA FILE REFERENCED WITH À 
P. 
ο πε”: ТИЕ RESUPTANT A PDLOT CANO PE DIRECIED TO ONE OF THREE 
PLOTTERS: VERSATEC 8222A, VERSATEC O900A, OR A 
TEKTRONIX 4631 HARDWARE DEVICE. 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


ооо о с СУ ОО ООО 


ШИШЕ ЗЕК Н А BTC DETA, AA IPETO ο 5 Т1МЕ,КА\МЕ,СОРЕ,АМТ,ВА5ЕТ 
1 SSA,SSB,SSC,SSD,SSE,ANTA,IPAK,ANTB,ANTC,ANTD,ANTE,SSX 


DIMENSION BASEL(1500), 

ANT(1500),CODE(1500),TIME(1500) Be ey 1500) 
SS(1500) , XRERR(1500) , XRATE(1500$, VRERR(1500) , VRATE(1500), 
XRERRX(1500) , SRERRA(1500) , XRERRB(1500 ) RRERRC(1500) , KRERRD(1500) , 
SSA(1500) ,SSB( 1500) ,SSC(1500) /SSD(1500) SSE (1 00) ,RAVE(1500), 
TEHPA(1500) TEMPB(1500) TEIPC 1500), TEMPD(1500),TEMPE(1500), 
SSX(1500),TEMPX(1 ΝῊ I 00500 ва 00), 
πο TOTALE (1500), S 1500),T(1500 
TOTALC(1500) , TOTALD(1500) , TOTALE(1500) 

REAL XRERR, XRERRA, XRERRB, XRERRC, XRERRD, XRERRE , XRATE, VRERR, 

1 TOTALA,TOTALB,TOTALC, TOTALD, TOTALE, TEMPX, 

1 TEMPA,TEMPB,TEHPC,TEHMPD, TEMPE , TOTAL, XRERRX 


3680р IN THE STAT DATA FROM A FIDEDEF. 


I-0 
7 I=I+1 
t 208) SS(I),XRATE(I),VRATE(I),XRERR(I),VRERR(I),BASEL(I), 
1СОВЕ(І) АМТ(І),БАУЕ(І),ТІМЕ(І) 
ІЕ(55(1).Е0.99% СО ТО 50 
Gamo. 7 


50 CONTINUE 
N-I 








, 


PB = 


(220) 


大 大 大 类 大 大 大 大 大 大 AT THIS POINT, ALL THE DATA HAS BEEN STORED IN ARRAYS. 
kkkkkkkkkk CALCULATE TOTAL MEAN FOR ALL DATA/CODE. 


71-0 
DO 31 M= 4,18 
TOTAL (M)=0 
DO 30 I=1,N 
πο. GO TO 30 
IF(SS(I).NE.M) GO TO 30 
TOTAL (M)=TOTAL(M)+XRERR(I) 


Х=Х+1 
30 CONTINUE 
XRERRX (M)=TOTAL(M) / FLOAT (X) 
SSX(M)=M 
ХХ=ХХ+Х 
31 CONTINUE 


C 
C ********^* DETERMINE, BY BASE LINE, THE MEAN RANGE ERRORS. 
С 


ОООО О 


ΠΩ 
H H H H l 
eOOOO 
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60 


70 


80 


90 


DO 165 M=4,18 
P (H 


Q (M 
R(M 
5(М 
T(H 
TOTALA | 


H Hl H| Il 
OOOOO 


TOTALB(M 
TOTALC(M 
TOTALD(M 
TOTALE (HM 
DO 160 I 


OOOOO 


CODE(I).NE.0) GO TO 160 
SS(I).NE.M) GO TO 160 
IF(BASEL(I).NE.1) GO TO 60 
Ai =TOTALA(M)+XRERR(I) 
P(M)=P +1 

GO TO 160 

IF(BASEL(I).NE.2) GO TO 70 
ОЕ E =TOTALB(M)+XRERR(I) 
о(и)= О М ii 


20 BASELII). NE.3) GO TO 80 
TOTALC(M) TOTALC(M) *XRERR(I) 
R(M)=R(M)+1 

GO TO 160 

IF(BASEL(I).NE.4) GO TO 90 
nw nc -TOTALD(M)-*XRERR(I) 
5(М)=5 +1 

GO TO 160 

ou S =TOTALE (M)+XRERR(I) 
T(M)=T(M) +1 


Mere jl i wou I 
"тј ету 
= 





160 CONTINUE 


SSA(M)TSS(I 
ο σος τ 
SSC(M)=SS(1 
SSD(M)-SS|I 
ӨЗЕТПЛЕЗ8 СІ 


C 
C 太 大 大 大 大 大 大 大 大 大 大 大 大 大 NOW CALCULATE THE MEAN RANGE ERROR FOR BAGH SIGNAL 


C 
С 


СуСУСУСУСУСУСУСУСУСУ 


O00 


STRENGTH FOR EACH BASE LINE. 


XRERRA (M)2TOTALA (M) /FLOAT(P(M 
XRERRB (M)-TOTALB(M)/FLOAT(Q(M 
XRERRC(M)=TOTALC(M)/FLOAT(R(M 
XRERRD(M)=TOTALD(M)/FLOAT(S(M 
XRERRE (M)=TOTALE(M) /FLOAT(T(M 
A=A+P (M 
B=B+0(M 
C=C+R(M 
D=D+S (M 
ЕзЕзТ(М 


165 CONTINUE 


大 大 大 大 大 大 大 大 大 大 大 大 大 PERFORM THE NECESSARY PLOT/SETUP ROUTINES. 
A***kkk**k INITIALIZE PLOTTER AND SET PAGE SIZE. 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


PRIPLT (72,6) 


COMPRS 
VRSTEC (0,0,0) 
TEK618 
PAGE (11.,8.5) 


BLOWUP(1. 55} 


大 大 火炎 大 炎炎 大 类 大 MAKE SS AXIS IN INTEGERS 
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anaa 


О0О О 


ONAA AAAA 


аса 


(NE 


CALL XINTAX 


в а ве оо аена -ᾱ 


ка 4%. ха 


SeGT ТОЛАТТРГОТ” SRBM 
CALL AREA2D(8.,6.) 


жилы С МА AAND LABEL AXES WITH SELP-COUNTING OPTIOM. 


CALL YNAME('RANGE ERROR (ii) Που 
CALL ANAME h SIGNAL STRENCIHS',100 
CALL BASALF('STAND' ) 


AXAXAAAA**** CREATE THE HEADING AND SET THE PRINT TYPE. 


ALL SWiSSM 
CALE HEEDIN ('MIMI-RANGER III BASE LINE GALIBRATIONS',100,1.5,1) 


° © a O στων „ә „ 3 


AAAAXA*X*** DEFINE BOTH AXES AND FRAME THE SUBPLOT AREA. 


CSEL DUPLX 
CALL CRAF(18.0, 
CALL THKFRM(.02 
CALI FRAM- 


1.03 „Ошаса тено в? 


CKDL LINESP (2.0) 

CALL LINES (!'CUMMULATIVES',IPAK,1) 
GREP SPSS (uk - 1061S‘, IPAK, 2 
οπως O DüNES (!35- 24175! ΙΡΆΚ,3 
СМ CINES ('С - 4083S' , IPAK,4 
САҒ. БІМЕ5 (!І) - 70175" ІРЕК /5 
CADL LINES ('B - 985601S', 3PAK,6 


9. 9. а орао еи ае. те 


ххх BIND STSE FORGDEGSND.. 


EE IAE 
YW-YLEGIIDD(IPAR,6 

CELL"MWAUBGN('BSASE LINE-LENGTH (M)S',100) 
CALL LEGLEN 


өлсек ы: ыла e 


KXX*X******* PLOT BASE LINE A VALUES. 


CALL MARKER(15) 
CALL DOT 
CALL TH IKCRV(. 025) 


Dee Ty's an 18 
TENeA(H) =F LOAT (SSA(M) ) 
195 CONTINUE 
CALL CUBNWE(TEMPA,XRERRA,R,1) 


Г р Т ЛАВОТ ТЕ е Т О δ- ë 


大 六 大 大 大 太太 大 六 六 六 PLOT BASE LINE B VALUES. 


CALL MARKER(15) 
CALL DASH 
CALL THKCRV(.025) 
DO 196 M=4,18 
TEMPB(M)=FLOAT(SSB(M) ) 
196 CONTINUE 
CALL CURVE(TEMPB,XRERRB,B,1) 


αι μὲ, ο. τ... να. ον че хиб са 0с ә а 


AXAXAAXA*XAAA** PLOT BASE LINE C VALUES. 


CABE MARRKER(I5) 
145 


ООО айя 


ООО О 


AQAA QAAAA ОООО) ООО 


ΩΩ 


CALL CHNDOT 
CALL THKCRV(. 025) 
DO 197 M=4,18 
TEMPC(M)=FLOAT(SSC(M) ) 
197 CONTINUE 
CALL CURVE(TEMPC,XRERRC,C,1) 


大 大 大 大 大 大 大 大 大 大 大 大 大 PLOT BASE LINE D VALUES. 


CALL MARKER(15) 
CALL CHNDSH 
CALL THKCRV(.025) 
DO 198 M=4,18 
TEMPD (M)=FLOAT(SSD(M)) 
198 CONTINUE 
CALL CURVE(TEMPD,XRERRD,D,1) 


WAkAkkkk*k*kk*k PLOT BASE LINE E VALUES. 


CALL MARKER(15) 
САП, МВ5СОМФМӘЕС 2.25. 2% 
CALL THKCRV(.025) 
DO 199 M=4,18 
TEMPE (M)=FLOAT(SSE(M) ) 
199 CONTINUE 
CALL CURVE(TEMPE , XRERRE,E,1) 


kkkkkkkkkkkk PLOT CUMMULATIVE VALUES. 
CALL MARKER(15) 
CALL RESET ('MRSCOD') 
CALL THKCRV(.025) 
DO 200 M=4,18 
TEMPX (M)=FLOAT(SSX(M) ) 
200 CONTINUE 
CALL CURVE(TEMPX,XRERRX,XX,1) 
kkkkkkkkkk OTHER PLOT/LEGEND INFORMATION: 


CALL SEEN TG МО, Ons Ιου. 
CALEMMESSAG('CODE NO. 15/7100 42550589 


大 大 大 大 大 大 大 大 大 大 DRAW THE LEGEND. 
CALL LEGEND ( Ρο ο ο 


大 大 大 大 大 大 大 大 大 大 READ FORMAT STATEMENT. 
208 FORMAT(6X,I2,7X,F6.1,8X,F3.1,10X,F4.1,11X,F3.1,3X,11,I11,11,12,14) 
大 大 大 大 大 大 大 大 大 大 CLOSE OUT METAFILE AND RETURN TO LEVEL ZERO. 

CALL ENDPL(0) 

CALL DONEPL 


SIOP 
END 
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